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ABSTRACT 


Because  of  the  complex  flow  fields  and  component  interaction  in  a  modem 
gas  turbine  engine,  these  engines  require  extensive  testing  to  validate  performance 
and  stability.  The  testing  process  can  become  expensive  and  complex.  Modeling 
and  simulation  of  gas  turbine  engines  is  one  way  to  reduce  testing  costs,  provide 
fidelity,  and  enhance  the  quality  of  essential  testing.  Several  numerical 
simulations  for  gas  turbine  engines  have  been  developed  at  Arnold  Engineering 
Development  Center  to  simulate  gas  turbine  engines  and  their  various  components. 
Compressor  performance  characteristics  are  needed  in  these  codes  to  provide 
turbomachinery  source  terms.  These  terms  are  currently  provided  by  experimental 
data.  A  technique  to  analytically  create  these  characteristics  would  greatly 
enhance  the  quality  and  value  of  the  codes  in  existence. 

Therefore,  a  technique  to  create  1-D  and  2-D  compressor  performance 
characteristics  with  specifications  of  only  annulus  geometry,  blade  geometry,  and 
loss  and  deviation  correlations  has  been  developed.  This  method  uses  a  previously 
developed  streamline  curvature  code  (SLCC)  with  open  literature  loss  and 
deviation  correlations  to  provide  on-  and  off-design  stage  performance.  Data 
reduction  techniques  are  then  used  to  convert  the  predicted  flow  field  behavior  in 
the  bladed  regions  into  stage-by-stage  performance  characteristics  of  the 


compressor. 


In  the  present  investigation,  it  was  discovered  that  the  SLCC  could  provide 
an  accurate  flow  field  prediction  with  the  calibration  of  relative  total  pressure  loss 
coefficients  and  flow  angle  deviations  to  experimental  data  and  the  inclusion  of 
inlet  and  exit  radial  blockages  to  correct  static  pressures  and  static  temperatures. 
The  results  of  the  1-D  representations  of  compressor  characteristics  for  a  calibrated 
flow  field  were  typically  within  one-percent  of  the  experimental  data.  The  2-D 
representation  gave  a  maximum  error  of  about  five-percent. 

It  was  also  discovered  that  machines  with  significantly  different  geometry 
will  have  different  trends  in  calibration.  A  second  machine  was  analyzed  with 
unfavorable  results  when  calibrated  with  the  trends  discovered  from  the  first 
machine.  This  second  machine  was  then  analyzed  without  calibrations  and 
produced  a  maximum  error  of  less  than  seven-percent  for  the  1-D  characteristics 
and  less  than  24-percent  for  a  2-D  representation. 
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1.0  INTRODUCTION 


Modem  aircraft  gas  turbine  engines  are  extremely  complex  machines, 
comprised  of  several  major  components.  One  of  the  major  components  of  the 
turbine  engine  is  the  compressor.  The  compressor  is  designed  to  increase  the 
pressure  of  the  incoming  air  and  deliver  the  correct  amount  of  airflow  to  the  rest  of 
the  engine.  A  common  type  of  compressor  for  large  gas  turbine  engines  is  the 
axial  flow  compressor.  An  axial  flow  compressor  consists  of  one  or  more  rotating 
blade  rows  which  are  aerodynamically  coupled  and  rotate  around  a  common  axis 
of  rotation.  A  typical  gas  turbine  compression  system  is  shown  in  Figure  l.l1.  A 
compressor  stage  consists  of  a  rotating  blade  row,  the  rotor,  and  a  stationary  blade 
row,  the  stator.  Rotor  blades  are  a  circumferential  array  of  airfoils  that  increase 
the  kinetic  energy  of  the  flow.  A  pressure  rise  across  the  blades  is  achieved  by 
flow  diffusion  in  the  rotor  blade  passage.  Stator  blades  diffuse  the  flow  more  and 
redirect  it  to  the  next  stage.  The  compressor  is  powered  by  a  combustor  and 
turbine  at  the  back  of  the  engine.  The  combustor  adds  energy  to  the  flow  and  the 
turbine  extracts  energy  from  the  flow  as  shaft  work  and  transfers  it  to  the 
compressor  by  a  shaft. 

Because  of  the  adverse  pressure  gradient  across  the  compressor,  its 
performance  and  stability,  on-  and  off-design,  are  critical.  The  performance  of  the 
compression  system  is  normally  shown  on  compressor  performance  maps.  These 

1  All  figures  and  tables  may  be  found  in  the  appendices. 
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maps  can  be  shown  in  the  form  of  pressure  ratio  and  efficiency  or  pressure  ratio 
and  temperature  ratio  versus  corrected  mass  flow.  Corrected  mass  flow  is  mass 
flow  which  has  been  adjusted  to  reference  total  pressure  and  total  temperature 
conditions.  Correcting  to  reference  conditions  collapses  the  performance  maps  for 
different  inlet  total  pressure  and  inlet  total  temperature  conditions  into  one  map. 
The  performance  data  presented  in  Figure  1.2  is  pressure  ratio  and  temperature 
ratio  versus  corrected  mass  flow.  The  maps  show  the  relationship  between 
compressor  total  pressure  ratio,  compressor  total  temperature  ratio,  corrected  mass 
flow,  and  corrected  engine  speed.  The  dashed  line  in  the  upper  part  of  the  total 
pressure  ratio  map  is  the  surge  line,  or  stall  line.  This  line  represents  the  limit  of 
the  pressure  ratio  as  a  function  of  corrected  rotor  speed  and  corrected  mass  flow 
rate.  To  the  right  of  the  line,  the  compressor  operates  in  a  stable  manner.  To  the 
left  of  the  stall  line,  the  compressor  may  not  operate,  or  will  operate  unstably. 
This  is  caused  by  high  enough  loading  on  the  blades  such  that  boundary-layer 
separation  occurs  over  a  large  part  of  the  blade.  Operation  beyond  the  stall  region 
is  structurally  dangerous  to  the  gas  turbine  engine  (Mattingly,  1996). 

Because  of  the  complex  flow  fields  and  component  interaction  in  gas 
turbine  engines,  these  engines  require  extensive  testing  to  validate  performance 
and  stability.  This  testing  process  is  very  expensive.  Furthermore,  it  is 
economically  and  technically  infeasible  or  impractical  to  measure  every  through- 
flow  quantity  at  every  location  throughout  the  engine.  Thus,  complete  component 
analysis  cannot  be  performed  during  testing.  Modeling  and  simulation  of  the  gas 
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turbine  engine  is  the  way  to  analyze  engine  performance  and  ensure  adequate  test 
fidelity  in  the  test  program. 

At  Arnold  Engineering  Development  Center  (AEDC),  modeling  and 
simulation  tools  are  looked  at  as  necessary  analytical  methods  to  reduce  testing 
cost  and  to  enhance  the  quality  of  essential  testing.  To  simulate  the  turbine  engine 
and  its  components,  several  gas  turbine  modeling  tools  have  been  developed  at 
AEDC.  These  codes  are: 

•  DYNamic  Turbine  Engine  Compressor  Code  -  DYNTECC 

•  Aerodynamic  Turbine  Engine  Code  -  ATEC 

•  Turbine  Engine  Analysis  Compressor  Code  -  TEACC. 

DYNTECC  (Hale  and  Davis,  1992)  and  ATEC  (Garrard,  1995)  are  one¬ 
dimensional  (1-D),  stage-by-stage  simulations.  DYNTECC  models  only  the 
compression  system  and  combustor,  while  ATEC  simulates  the  entire  gas 
generator  (compressor,  combustor  and  turbine).  The  governing  equations  in  both 
simulation  models  are  the  mass,  momentum,  and  energy  conservation  equations 
applied  to  elemental  control  volumes  with  turbomachinery  source  terms  (mass 
bleed,  blade  forces,  heat  transfer,  and  shaft  work)  included  in  the  blade  row 
elemental  control  volumes.  The  source  terms  simulate  the  effect  the  blades  have 
on  a  control  volume  and  are  determined  from  a  complete  set  of  stage  pressure  and 
temperature  characteristics. 
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TEACC  (Hale  and  O’Brien,  1997)  is  a  three-dimensional  (3-D)  compressor 
code.  This  code  solves  the  compressible,  3-D  Euler  equations  modified  to  include 
turbomachinery  source  terms.  Again,  compressor  characteristics  can  be  used  to 
calculate  these  source  terms. 

Currently,  the  performance  characteristics  must  be  provided  by 
experimental  data  for  a  given  stage  or  compressor.  To  be  able  to  model  a  machine 
without  experimental  data,  a  general  method  is  needed  to  generate  the  performance 
maps  for  a  given  compressor.  The  work  included  in  this  thesis  will  provide  the 
necessary  steps  of  using  compressor  geometry  parameters,  a  streamline  curvature 
code  (SLCC),  and  data  reduction  to  develop  both  1-D  and  radial  distributions  of 
performance  characteristics  to  calculate  the  turbomachinery  source  terms  required 
by  each  of  the  three  codes  mentioned  above. 

This  thesis  is  organized  to  first  give  a  brief  review  of  other  techniques  used 
to  develop  compressor  performance  characteristics.  The  methodology  and  solution 
technique  of  the  SLCC  are  then  explained.  This  section  includes  an  in-depth 
investigation  of  the  transonic  loss  and  deviation  correlations  used  in  the  SLCC. 
Next,  the  calibration  of  the  SLCC  to  a  machine  with  experimental  data  is 
presented.  With  the  results  of  the  calibration  for  the  first  machine,  the  1-D  and 
two-dimensional  (2-D)  performance  characteristics  of  another  machine  are 
predicted.  Then,  the  results  are  summarized  and  conclusions  are  discussed. 
Finally,  suggestions  for  future  work  are  presented. 
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2.0  REVIEW  OF  THE  LITERATURE 

Because  stage-by-stage  analysis  of  a  gas  turbine  compression  system 
requires  compressor  stage  performance  characteristics,  many  attempts  have  been 
made  to  develop  these  characteristics.  These  methods  range  from  a  1-D,  2-D,  and 
3-D  analyses,  to  altering  generic  maps  to  match  limited  stage-by-stage  data.  Prior 
to  explaining  the  current  methodology,  previous  work  on  generation  of 
performance  characteristics  is  reviewed. 

2.1  Experimental  Data 

Experimental  methods  are  a  very  common  way  to  develop  compressor 
performance  characteristics.  Pressure  and  temperature  measurements  are  obtained 
at  the  front  and  back  of  a  compressor  stage  (Cyrus,  1996).  This  data  can  be  taken 
from  either  compressor  rig  tests  or  actual  complete  engine  tests.  The  pressure  and 
temperature  data  is  then  reduced  to  stage-by-stage  characteristics.  These  methods 
require  the  actual  testing  of  hardware.  Also,  in  the  development  phase,  to  test  a 
new  blade  shape  or  stage  setup,  new  blades  have  to  be  cut  and  tested.  Davis  and 
O’Brien  (1991)  took  experimental  data  from  a  three-stage  compressor  rig  and 
converted  it  to  1-D,  steady-state  characteristics.  These  steady-state  characteristics 
were  then  used  to  model  an  unsteady,  post-stall  compression  system.  Similar 
studies  were  presented  by  Davis,  Hale,  Sharohroki,  and  Garrard  (1996)  for  a  10- 
stage  machine  and  a  T-55  helicopter  engine.  Because  these  methods  require  the 
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building  and  testing  of  hardware,  they  can  become  expensive.  Therefore  accurate 
computational  methods  are  desirable. 

2.2  One-Dimensional  Techniques 

One-dimensional  (1-D)  techniques  are  simple  and  fast  methods  for 
obtaining  compressor  performance  predictions.  A  drawback  to  the  1-D  analysis  is 
that  the  true  3-D  physics  of  the  flow  in  the  compressor  are  replaced  by  a  1-D 
model  which  must  inherently  predict  the  average  of  what  occurs  across  a  blade 
row.  A  1-D  method  to  generate  stage  characteristics  was  presented  by  Attia  and 
Schobeiri  (1995).  This  method  required  that  blade  row  exit  flow  angles  be  given 
and  used  a  loss  correlation  model  developed  by  one  of  the  authors  to  predict  loss 
across  a  blade  row.  The  only  other  required  inputs  were  geometry,  rotational 
speed,  and  mass  flow  rate.  The  technique  assumed  that  deviation  angles  were 
negligible.  The  deviation  angle  is  the  angular  difference  between  the  blade  metal 
angle  and  the  actual  exit  flow  angle.  The  results  were  presented  in  the  form  of 
pressure  ratio  and  efficiency  versus  mass  flow.  It  was  concluded  that  the  model 
was  satisfactory  for  blades  of  small  height.  This  result  were  satisfactory  only  for 
blades  of  small  height  because  of  the  assumption  that  the  deviation  angle  was 
negligible  everywhere,  which  cancelled  hub  and  tip  effects. 

Johnson  (1991)  presented  a  method  for  calculating  1-D,  stage-by-stage 
characteristics  using  a  blade  element  technique.  This  technique  requires  the  actual 
inlet  and  exit  flow  conditions  and  solves  the  continuity,  energy,  and  momentum 
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equations  at  every  stage.  He  used  a  deviation  correlation  presented  by  Horlock 
(1958).  Horlock  (1958)  also  provided  a  relationship  between  coefficient  of  drag 
and  total  pressure  loss  coefficient  that  was  used  in  this  correlation.  The  data  used 
for  both  the  loss  and  deviation  correlations  were  developed  from  NACA  data 
collected  by  Emery  (1958).  For  the  compressor  studied  by  Johnson,  this  technique 
yielded  predictions  for  performance  characteristics  that  were  higher  than 
experimental  data.  This  inaccuracy  was  blamed  on  the  total  pressure  loss 
coefficient  correlation. 

The  final  1-D  method  cited  here  was  by  Tsalavoutas  et.  al.  (1994).  This 
method  was  essentially  a  stage  stacking  technique.  It  required  no  knowledge  of 
internal  geometry  but  required  some  experimental  performance  data  for  the 
compressor.  This  method  also  involved  no  flow  physics  and  relied  entirely  on 
experimental  data.  The  method  started  with  an  initial  set  of  generic  stage 
characteristics  and  used  optimization  routines  until  selected  points  on  the  overall 
map  were  obtained  within  a  desired  accuracy  level.  The  individual  stage 
characteristics  may  not  have  been  accurate,  but  the  overall  map  was.  The  results 
of  the  characteristics  were  in  the  form  of  pressure  ratio  and  efficiency  versus 
corrected  mass  flow. 

2.3  Two-Dimensional  Techniques 

Two-dimensional  techniques,,  or  throughflow  techniques,  are  another 
common  way  to  develop  characteristics.  Sayari  and  Boies  (1995)  developed  a 
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technique  based  on  the  streamline  curvature  method.  Most  streamline  curvature 
codes  (SLCC)  use  a  circumferential  averaging  technique  to  develop  the  2-D  flow 
field.  Sayari  and  Boies  offered  a  method  that  averaged  normal  to  the  central 
streamline.  They  argued  that  this  method  more  accurately  predicts  the  Mach 
number  just  before  and  around  the  shock.  This  model  was  considered  an 
improvement  on  the  commonly  used  Miller  shock  model  discussed  in  detail  in 
Section  3.2.1.  This  method  concentrated  on  loss  and  deviation  due  to  shocks  and 
gave  no  consideration  to  other  loss  and  deviation  possibilities.  The  machines 
predicted  were  only  simulated  in  the  region  where  shocks  were  present.  They 
presented  graphical  comparisons  to  data  but  never  mentioned  percent  differences. 

Korakianitis  and  Zou  (1993)  presented  a  streamline  curvature  method  that, 
instead  of  iteratively  solving  the  radial  momentum  equation,  used  a  predictor- 
corrector  technique  to  solve  it  in  one  pass.  This  reduced  computational  time.  No 
correlations  for  deviation  were  presented;  they  must  be  provided  by  the  user  as 
known  inputs.  Loss  across  a  blade  row  could  have  either  been  a  user  input  or 
calculated  by  using  open-literature  loss  correlations. 

The  final  2-D  method  (Yazigi,  N.  et.  al.,  1990)  reviewed  used  a  streamline 
curvature  method  for  performance  prediction.  This  method  incorporated  an 
inverse  boundary-layer  model  to  calculate  the  boundary-layer  thickness  at  the 
trailing  edge.  Scholz’s  model,  as  described  by  Klein  (1977),  was  used  for 
estimation  of  the  mixing  losses.  Deviation  was  determined  by  potential  flow 
theory,  taking  into  account  the  influence  of  the  boundary  layers,  separation,  and 
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stream  tube  height.  The  performance  was  presented  in  the  form  of  pressure 
distribution.  The  computational  time  required  for  this  method  was  still  quite  high, 
at  nearly  0.75-hours  on  an  IBM-PC  for  one  run. 

2.4  Three-Dimensional  Techniques 

Another  computational  method  of  developing  performance  characteristics 
is  the  use  of  3-D  computational  fluid  dynamics  (CFD)  codes.  A  method  presented 
by  Hah  and  Wennerstrom  (1990)  used  a  Reynolds-averaged  Navier-Stokes 
equation  solution  to  develop  the  flow  field  in  blade  rows.  The  results  of  the  3-D 
flow  field  for  a  transonic  compressor  with  swept  blades  were  presented  as  total 
pressure  and  total  temperature  characteristics.  The  resulting  characteristics  were 
within  2-percent  of  experimental  data.  This  code  required  one  hour  of 
computational  time  on  a  Cray  XMP  supercomputer  to  converge  to  a  solution  for 
one  flow  point. 

Another  three  dimensional  technique  was  presented  by  Swan  (1961).  This 
technique  used  non-isentropic  complete  radial  equilibrium-momentum  equations 
combined  with  continuity,  energy,  and  flow-process  relations  to  derive  a  group  of 
equations  that  could  be  solved  on  a  computer.  The  results  of  a  single  stage 
machine  were  presented  in  the  form  of  total  pressure  ratio  and  adiabatic  efficiency 
versus  corrected  mass  flow.  The  results  given  were  accurate  to  within  10-percent 
of  experimental  data.  This  method  used  only  a  single  variable  to  statistically 
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deduce  the  real  fluid  effects  such  as  viscous  losses  and  shock  losses.  The  method 
was  computer  time  intensive. 

2.5  Present  Investigation 

As  a  consequence  of  the  brief  review  of  these  previous  compressor 
modeling  approaches,  it  was  concluded  that  the  previous  works  lacked  important 
desirable  features  need  for  this  study.  Specifically,  all  models  exhibited  one  or 
more  of  the  following  deficiencies: 

•  Experimental  data  for  the  particular  machine  under  study  is  required 

•  Only  one  or  two  sources  of  loss  and  deviation  are  utilized  in  the  model 

•  Model  does  not  account  for  spanwise  variation  in  the  flow 

•  Computational  times  are  too  long 

The  purpose  of  this  thesis  is  to  develop  a  compressor  performance 
prediction  technique  that  overcomes  the  deficiencies  of  the  previous  models  and 
can  give  sufficiently  accurate  compressor  performance  predictions  for  compressors 
operating  in  pre-stall  conditions.  The  technique  will  be  based  on  input  of 
geometric  data  for  the  annulus  and  blade  shape;  no  experimental  data  for  the 
particular  machine  under  study  is  needed,  except  for  calibration  and  validation.  It 
is  based  on  an  existing  streamline  curvature  code  with  several  general  loss  and 
deviation  correlation  models.  The  model  gives  a  radial  distribution  of  the  entire 
flow  field  and  requires  only  seconds  to  run  a  flow  point  on  a  desktop  computer. 
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3.0  APPROACH 


The  approach  used  in  the  present  investigation  is  based  on  the  use  of  a 
streamline  curvature  code  (SLCC).  The  SLCC  is  essentially  a  duct  flow  solver 
with  an  imbedded  set  of  loss  and  deviation  correlations  which  provide  the  means 
to  step  across  a  bladed  region.  Because  the  SLCC  is  an  integral  part  of  this 
approach  to  calculating  compressor  stage  characteristics,  and  the  correlations  are 
the  heart  of  the  SLCC,  each  needs  to  be  investigated.  A  data  reduction  technique 
was  developed  that  converts  a  radial  distribution  of  performance  quantities  to  their 
equivalent  1-D  representations,  and  generates  performance  characteristics  from  the 
SLCC  output.  Details  of  the  data  reduction  technique  are  also  presented  in  this 
thesis. 

3.1  Streamline  Curvature  Code  Methodology 

In  the  early  1950’s,  Wu  (1952)  developed  the  governing  equations  and 
general  theory  for  a  three-dimensional  analysis  of  turbomachines.  With 
improvements  in  computer  sizes  and  speeds  in  the  mid-1960’s,  numerical  solution 
methods  based  on  Wu’s  general  theory  were  developed  to  solve  flows  through 
turbomachines.  Marsh  (1966)  developed  the  matrix  through-flow  method,  and 
Smith  (1966)  developed  the  streamline  curvature  method.  Streamlines  are  defined 
as  lines  which  are  tangent  to  the  velocity  vectors  throughout  the  flow  field.  The 
method  used  in  the  current  approach  to  determine  the  flow  field  is  a  streamline 
curvature  method  developed  by  Hearsey  (1994)  and  is  similar  to  the  one  developed 
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by  Smith  (1966).  The  SLCC  is  a  computer  program  that  solves  an  axisymmetric, 
annular  flow  field.  The  flow  is  assumed  to  be  an  inviscid,  perfect  gas  with  no 
transfer  of  mass,  momentum,  or  energy  between  adjacent  streamlines. 

The  inputs  required  by  the  SLCC  are  the  overall  annulus  geometry,  blade 
geometry,  and  a  radial  distribution  of  inlet  total  temperature,  total  pressure,  swirl 
angle,  and  loss  and  deviation  correlations  which  are  sensitive  to  the  local  flow 
field.  The  annulus  geometry  is  divided  into  a  series  of  axial  stations  that  need  not 
necessarily  be  radial  but  may  be  leaned,  or  inclined,  at  an  angle  with  respect  to  the 
radial  direction  to  facilitate  placement  at  locations  of  interest,  such  as  the  leading- 
and  trailing-edge  blade  angles.  Initial  estimated  streamlines  are  then  set  up 
radially  to  create  a  computational  grid.  The  slope  of  the  streamlines  are  assumed 
to  be  zero  far  from  the  bladed  regions  at  the  inlet  and  exit.  A  schematic  of  the 
SLCC  inputs  and  how  the  annulus  is  divided  into  a  computational  grid  is 
illustrated  in  Figure  3.1. 

A  computing  station  is  defined  by  the  computational  grid.  A  sketch  of  the 
computing  station  is  shown  in  Figure  3.2.  In  this  figure,  r  is  the  radial  direction,  z 
is  the  axial  direction,  /  is  the  computing  station  direction,  m  is  the  meridional 
direction,  rc  is  the  radius  of  curvature  of  the  streamline,  y  is  the  station  lean  angle, 
and  <pis  the  streamline  slope  angle. 

The  velocity  triangle  nomenclature  used  in  the  SLCC  is  presented  in  Figure 
3.3.  The  meridional  velocity  (Vm)  is  the  base  velocity,  defined  as  the  vector  sum 
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of  the  axial  and  radial  velocity  components.  One  limitation  of  the  streamline 
curvature  method  (Smith  1966)  is  that  the  velocities  in  the  absolute  reference  plane 
must  be  subsonic.  The  relative  Mach  numbers  have  no  restrictions.  Alpha  (a)  is 
the  absolute  flow  angle  and  is  defined  as  the  angle  between  Vm  and  the  absolute 
velocity,  V.  Beta  (ft)  is  the  relative  flow  angle  and  is  defined  as  the  angle  between 
Vm  and  the  relative  velocity  (W).  U  represents  the  wheel  rotational  velocity. 

The  following  paragraphs  briefly  describe  the  solution  technique  of  the 
SLCC.  A  complete  derivation  of  the  governing  equations  can  be  found  in 
Appendix  C.  The  SLCC  assumes  an  inviscid,  axisymmetric,  adiabatic,  steady- 
state  flow  with  no  body  forces.  The  governing  equations  are  continuity, 
circumferential  momentum,  axial  momentum,  radial  momentum,  energy,  entropy 
and  the  ideal  gas  equation  of  state.  Mapping  the  circumferential  momentum 
equation  to  the  non-orthogonal,  m  and  l  coordinate  system  reveals  that  angular 
momentum  is  constant  along  a  streamline.  As  shown  in  Appendix  C,  mapping  of 
the  radial  and  axial  momentum  equations  to  the  m  and  /  coordinate  system,  in 
combination  with  the  inclusion  of  the  T-dS,  energy,  and  entropy  equations, 
reveals  that  the  radial  and  axial  momentum  equations  are  equivalent.  For  this 
reason,  the  axial  momentum  equation  is  discarded. 

Duct  flow  regions  and  bladed  regions  are  solved  by  different  equations. 
Both  methods  use  the  radial  momentum  equation,  continuity  equation,  ideal  gas 
equation  of  state,  definition  of  static  temperature,  and  the  isentropic  relationship 
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for  static  pressure.  The  other  equations  involved  are  used  to  calculate  three 
closure  parameters.  These  parameters  are  exit  swirl  velocity,  exit  total 
temperature,  and  exit  total  pressure. 

The  closure  relations  for  duct  flow  regions  are  presented  first.  From 
integration  of  the  circumferential  momentum  equation,  energy  equation,  and 
entropy  equation  along  a  streamline,  it  is  revealed  that  angular  momentum, 
enthalpy,  and  entropy  are  constant  along  a  streamline.  This  fact  gives  the  three 
closure  relations  as 


Vn  = 


vr* 


r  ev 


Tn  —  , 

^12  =  Pfl  ’ 


[3.1] 

[3.2] 

[3.3] 


where  the  subscripts  1  and  2  represent  the  inlet  and  exit  to  a  computational  station. 

For  a  bladed  region,  the  process  is  more  complex.  Angular  momentum, 
enthalpy,  and  entropy  are  not  constant  along  a  streamline.  Because  of  the 
complexity  of  the  flows  involved  in  the  bladed  region,  the  SLCC  does  not  directly 
model  the  bladed  region.  Empirically  derived  loss  and  deviation  correlations 
provide  the  relative  total  pressure  loss  coefficient  and  flow  deviation  across  a 
bladed  region.  The  relative  total  pressure  loss  coefficient  is  defined  as  the  mass- 
averaged  defect  in  relative  total  pressure  divided  by  the  pressure  equivalent  of  the 
inlet  velocity  head, 
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[3.4] 


i=  Vt2  )jd  ^t2 

P'n-P, 

where  (pn  )jd  is  the  ideal  exit  relative  total  pressure  with  no  losses,  Pt2  is  the  exit 
relative  total  pressure  with  losses,  and  the  denominator  is  the  relative  inlet 
dynamic  pressure  head. 

The  deviation  correlations  provide  the  deviation  of  the  exit  flow  angle,  fo, 
from  the  blade  metal  angle,  /%m, 

<?  =  &-&„•  P-5] 


The  correlations  for  both  loss  and  deviation  are  experimentally  derived  and 
described  in  detail  in  Section  3.2. 

With  loss  and  deviation  provided  by  the  correlations  at  each  streamline, 
the  closure  relations  that  specify  the  change  in  enthalpy,  entropy,  and  angular 
momentum  across  a  blade  row  are  developed.  The  relations  are  derived  in 
Appendix  C.  The  final  relations  are 


vd2=vm2  tan  f32  cos +U2, 
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With  the  exit  swirl  velocity,  exit  total  temperature,  and  exit  total  pressure 
determined,  the  conservation  equations  described  in  Appendix  C  for  each  case  are 
solved  to  develop  the  entire  flow  field. 

3.2  Correlations  Used  in  the  Streamline  Curvature  Code 

Because  of  the  complex  physical  flow  occurring  in  a  compressor  stage,  it  is 
very  difficult  to  model  the  flow  phenomena  exactly.  For  this  reason,  empirical 
correlations  are  typically  used  to  approximate  the  actual  physics  involved.  Most  of 
the  correlations  used  in  the  SLCC  were  originally  derived  from  2-D  linear  cascade 
flow  results  and  are  found  in  NASA  SP-36  (1965).  These  correlations  were 
developed  from  databases  acquired  from  machines  of  1950’s  and  1960’s  design. 
Because  they  were  developed  on  earlier  designs,  the  correlations  may  have  trouble 
when  modeling  machines  with  more  modem  blade  profiles.  However,  some 
correlations  have  been  modified  to  better  approximate  modem  high  speed 
turbomachines. 

Before  the  loss  and  deviation  correlations  can  be  used,  several  geometric 
quantities  must  be  known.  These  are  the  blade  inlet  and  exit  metal  angles,  solidity, 
camber,  blade  maximum  thickness  to  chord  ratio,  and  location  of  maximum 
camber  point  as  a  fraction  of  chord.  With  these  quantities  defined,  two  reference 
values  are  calculated  using  the  equations  found  NASA  SP-36  (Lieblein,  1965). 
The  first  parameter  is  the  low-speed  reference  minimum-loss  incidence  angle 
given  by 
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Kef  =  (ki  ) shape  (*i  \h  (*0  )lO  “ nd  ’  [3'9] 

where  (k/)ShaPe  and  (k,)th  are  correction  factors  for  blades  with  shapes  other  than 
NACA  65-(Aio)-series  blades  and  (i0)io  is  the  variation  of  zero-camber  incidence 
angle  for  the  10-percent-thick  65-series  thickness  distribution,  (r’o)io  is  a  function 
of  inlet  air  angle,  fi\.  This  value  is  unknown  at  the  time  of  calculation  of  iie f,  so  it 
is  found  iteratively.  The  term  n  is  the  minimum-loss-incidence  slope  factor.  irefis 
used  in  Equations  [3.10]  and  [3.34]  -  [3.37]. 

The  last  reference  value  that  must  be  found  is  the  reference  minimum-loss 
air  inlet  angle.  Assuming  the  minimum-loss  incidence  angle  is  zero,  the  reference 
minimum-loss  air  inlet  angle  would  be  equal  to  the  inlet  blade  metal  angle. 
However,  Equation  [3.9]  calculates  a  reference  minimum-loss  incidence  angle,  so 
a  better  initial  estimate  of  the  reference  minimum-loss  air  inlet  angle  is 

/W-A.+'*-  P-io] 

p\ref  is  used  in  Equations  [3.13],  [3.15],  and  [3.23]. 

3.2.1  Loss  Correlations 

The  loss  correlations  calculate  the  relative  total  pressure  loss  coefficient, 

W' ,  defined  by  Equation  [3.4].  The  SLCC  calculates  loss  at  a  given  compressor 
operating  point  by  developing  a  "loss  bucket".  A  loss  bucket  is  the  graphical 
representation  of  loss  as  a  function  of  incidence  for  constant  Mach  number.  The 
curve  normally  has  a  U-shape  or  bucket  shape  with  the  middle  of  the  bottom  of  the 
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bucket  being  the  minimum-loss  incidence  point.  The  equation  that  develops  the 
loss  bucket  for  the  SLCC  is 

®'=  (®«nn+®M+®Hub+®Jl  +  (1-1^)  ■  [3.11] 

L  v  w  J  J. 

The  first  four  terms  are  magnitudes  of  different  loss  values,  i  is  the  actual 
incidence,  i'm  is  the  minimum  loss  incidence,  and  W  is  an  arbitrarily  defined  width 
of  the  loss  bucket.  iM  is  a  function  of  inlet  relative  Mach  number  and  iK f  (defined 
by  Equation  [3.9]).  An  example  of  the  loss  bucket  and  how  each  term  effects  the 
bucket  shape  is  shown  in  Figure  3.4.  The  solid  line  represents  a  baseline  bucket. 
The  baseline  bucket  is  an  arbitrary  bucket  that  will  be  used  as  a  baseline  to  depict 
how  the  other  variables  affect  the  bucket.  With  increasing  minimum  incidence 
(z'm),  the  bucket  shifts  to  the  right.  With  increasing  magnitude  of  loss  terms  (GW 
G5m,  GJhub,  GJtip),  the  bucket  shifts  up.  As  the  width  term  (W)  decreases,  the  bucket 
width  decreases. 

The  minimum  loss  term,  G7min ,  is  essentially  the  profile  loss,  or  friction  loss, 
caused  by  viscous  forces  on  the  blade.  The  one  used  for  this  investigation  is 
proposed  by  Robbins  et.  al.  (1965)  in  NASA  SP-36.  The  correlation  for  minimum 
profile  loss  is  given  as 


where  GJvloss  is  an  input  loss  modification  term.  The  P  term  in  Equation  [3.12]  is 

interpolated  from  a  data  set  developed  from  Figure  3.5.  The  abscissa  for  Figure 
3.5  is  the  design  diffusion  factor,  Ddes-  The  diffusion  factor  used  in  the  SLCC  is 
based  on  flow-weighted  mean  inlet  and  exit  axial  velocities  and  the  outlet  flow 
angle  plus  the  deviation  due  to  mean  axial  velocity  changes.  Equation  [3.13] 
shows  this  relationship 


2 (r,Vol  tan  filref  -  r2Va2  tm(fi2ref  +  Svamn ))  2 nN s 
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(r2~ri)  X 


2(0.45  +  0.5fe)_  )]  +  5.0fcL,  •  [3-13] 

2V.,C  J 

where  Va  is  the  axial  velocity,  &amn  is  the  deviation  due  to  mean  axial  velocity 
defined  similar  to  Equation  [3.21]  except  using  Va  instead  of  Vm,  and  Ns  is  the 
rotational  wheel  speed  in  terms  of  revolutions  per  minute  (rpm). 

The  next  loss  correlation  to  be  examined  is  the  loss  due  to  a  strong  shock  in 
the  passage.  The  model  for  shock  loss  is  the  one  presented  by  Miller,  Lewis,  and 
Hartman  (1960).  The  inlet  critical  Mach  number  is  defined  as  the  value  of  the 
inlet  Mach  number  that  will  give  supersonic  flow  somewhere  in  the  blade  passage. 
This  value  is  found  by  assuming  that  the  pressure  coefficient  corresponding  to  the 
minimum  pressure  point  on  the  blade  suction  surface  remains  practically 
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unchanged  up  to  the  critical  Mach  number  (Grewe,  1957).  From  this  assumption, 
the  following  equation  for  critical  inlet  Mach  number  was  derived  (Al-Daini, 


1986) 

[3.14] 


Because  Equation  [3.14]  requires  the  ratio  of  maximum  relative  velocity  to  inlet 
relative  velocity,  the  following  correlation  (Jansen  and  Moffat,  1967)  was 
developed 
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A  schematic  of  the  shock  loss  model  is  shown  in  Figure  3.6.  If  the  inlet 
relative  Mach  number  is  larger  than  its  critical  value,  a  shock  will  exist  somewhere 
in  the  blade  passage.  The  shock  will  be  a  detached  bow  shock  but,  to  simplify  the 
problem,  it  is  assumed  to  be  a  normal  shock.  The  normal  shock  extends  from  the 
tip  of  the  blade,  normal  to  the  mid-channel  streamline,  and  intersects  the  following 
blade  on  the  suction  side,  shown  in  schematically  in  Figure  3.6.  The  Mach  number 
at  point  A  in  Figure  3.6  is  assumed  to  be  the  inlet  relative  Mach  number.  The 
Mach  number  at  B  is  found  by  using  a  Prandtl-Meyer  expansion  from  the  blade  tip 
to  the  point  of  shock  impingement.  The  expansion  angle  is  approximated  by 
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assuming  the  blade  is  a  double-circular-arc  (DCA)  section.  The  shock  Mach 
number,  Ms,  is  an  average  of  the  inlet  relative  Mach  number  and  the  Mach  number 
at  point  B.  Solidity  is  important  in  this  correlation  because  blade  spacing  will 
directly  affect  the  point  of  shock  impingement,  hence  affect  the  Prandtl-Meyer 
expansion  angle.  With  the  shock  Mach  number  known,  the  loss  coefficient 
associated  with  the  shock  given  by  Lewis,  Miller,  and  Hartman  is  developed  from 
normal  shock  and  isentropic  relations  and  is  given  as 

i 

Y- 1 

-1 

- .  [3.16] 


The  last  two  loss  terms  in  Equation  [3.11]  are  the  losses  due  to  hub  and  tip 
effects,  such  as  intense  and  concentrated  vorticity  in  the  flow.  Because  the  rotors 
centrifuge  some  boundary-layer/wake  flow  radially  outwards,  increased  losses 
appear  at  the  tip  (Hearsey,  1994).  Hearsey  also  states  that  some  experimental  data 
shows  an  increase  in  the  loss  coefficient  near  the  hub  that  is  not  accounted  for  by 
the  other  correlations.  Therefore,  he  includes  the  0huh  term.  Both  GJtip  and  G7hub 

are  simply  cubic  functions  of  radius  that  add  loss  at  both  the  hub  and  tip  in 
proportion  to  user  inputs,  Hioss  and  Tioss-  The  equations  for  hub  and  tip  loss  are 
given  as 
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Hearsey  suggests  both  Hloss  and  Tioss  to  be  zero  for  stationary  blade  rows.  For  a 
rotor,  he  suggests  Hioss  be  set  to  match  experimental  data  and  Tioss  be  set  to  one  for 
the  first  rotating  blade  row  and  decreased  by  0.2  for  each  rotating  blade  row 
afterwards  until  it  also  is  zero. 

3.2.2  Deviation  Correlations 

The  deviation  correlations  represent  the  deviation  in  flow  angle  from  the 
actual  blade  trailing  edge  or  exit  metal  angle  (Figure  3.7).  Such  deviations  are 
caused  by  thick  or  separating  boundary-layers,  usually  on  the  suction  side  of  the 
blade.  Like  loss,  there  are  several  different  reasons  for  deviation  to  occur.  Five 
specific  reasons  will  be  investigated  in  this  section.  The  total  deviation  is  assumed 
to  be  representable  as  a  linear  combination  of  the  five  reasons.  The  total  deviation 


is  given  as 


8  =  $  ref  +  <>va  +  ^3  D  +  +  <>i  • 
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The  first  deviation  correlation  investigated  (Lieblein,  1965)  is  the  low- 


speed  reference  minimum  loss  deviation  angle 
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where  (ks)ShaPe  and  (kg)*  are  correction  factors  for  blades  with  shapes  other  than 


NACA  65-(Aio)-series  blades  and  (<3b)io  is  the  variation  of  zero-camber  deviation 
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angle  for  the  10-percent-thick  65-series  thickness  distribution.  (&)io  is  a  function 
of  inlet  air  angle,  / \ .  The  term  m  is  the  deviation  slope  factor  and  b  is  the  solidity 
exponent.  The  equation  for  varies  from  the  original  one  presented  by  Lieblein. 
A  modification  to  it  based  on  previous  empirical  data  (Carter,  1946)  made  cW  a 
function  of  (a/c),  the  point  of  maximum  camber  as  a  fraction  of  chord. 

The  next  deviation  correlation  is  based  on  axial  velocity  ratio  (Horlock, 
1967-68).  Horlock  states  that  because  of  the  growth  of  end  wall  boundary  layers, 
the  pressure  increase  (by  diffusion)  in  the  blade  row  is  decreased,  resulting  in  an 
axial  velocity  ratio  that  is  different  from  the  design.  Thus,  the  axial  velocity  ratio 
will,  in  general,  not  be  unity.  Figure  3.8  shows  how  the  streamline  contraction 
across  a  blade  row  causes  the  streamtube  area  to  change.  With  mass  not  allowed 
to  cross  a  streamline,  the  streamtube  area  contraction  requires  that  the  axial 
velocity  must  increase  (i.e.  diffusion  is  reduced).  With  the  axial  velocity  ratio 
larger  than  one,  he  noted  the  deviation  actually  decreases.  Therefore,  he  presented 
this  correlation 


The  third  deviation  correlation  is  the  deviation  due  to  streamline  radial 
location.  Because  of  complex,  three-dimensional  flows  at  the  blade  hub  and  tip, 
the  deviation  at  these  locations  can  be  different  from  the  deviation  at  the  mid-span. 
This  correlation  (Robins,  Jackson,  and  Lieblein,  1965),  <%D,  accounts  for  some  of 


23 


these  effects.  For  high-speed  sections,  or  blades  with  a  DCA  profile,  S3D  is 
interpolated  from  Figure  3.9.  In  this  figure,  S3D  is  a  function  of  the  inlet  relative 
Mach  number  and  percent  of  blade  height  from  compressor  tip.  For  low-speed 
sections,  or  blades  with  NACA  65-(A]0)-series  sections,  S3D  is  equal  to  -0.5 


degrees  for  blade  height  greater  than  50-percent  of  blade  height  from  compressor 
hub.  For  all  other  radial  locations  of  NACA  65-(Ai0)-series  sections, 


$3D  =“K+9% 


[3.22] 


The  existence  of  a  shock  in  the  blade  passage  is  another  source  of 
deviation.  This  correlation  is  only  used  if  the  inlet  relative  Mach  number  is  larger 
than  one.  The  development  of  this  correlation  involves  knowledge  of  the  inlet  and 
exit  critical  Mach  numbers  and  the  inlet  and  exit  axial  velocities.  The  inlet  critical 


Mach  number  is  presented  in  Equation  [3.14].  In  the  same  paper  (Jansen  and 
Moffat,  1967),  the  exit  critical  Mach  number,  the  inlet  and  exit  critical  velocities, 
and  the  inlet  and  exit  axial  velocities  are  derived.  The  exit  critical  Mach  number, 
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is  solved  iteratively  across  the  blade  and  used  to  calculate  the  exit  critical  velocity. 
The  exit  critical  Mach  number  is  also  used  to  calculate  the  exit  axial  velocity.  The 


inlet  and  exit  critical  velocities  given  by  Jansen  and  Moffat  are 
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The  actual  outlet  flow  angle  caused  by  the  inlet  Mach  number  above  the 


critical  value  is  given  by 


cos  fl2M  =  cos  /?, 
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The  deviation  due  to  the  Mach  number  greater  than  one  in  the  blade  passage  is 
then  simply 

P-3®] 

The  final  deviation  component  is  the  deviation  due  to  off-design  incidence. 
To  determine  this  value,  the  stalling  and  choking  incidences  must  be  known.  The 
stalling  incidence  is  the  incidence  at  which  a  boundary-layer  separation  occurs 
over  a  large  part  of  the  blade,  causing  the  blade  to  stall.  The  choking  incidence  is 
the  flow  incidence  at  which  the  flow  in  the  blade  passage  becomes  choked.  The 
equations  for  stalling  and  choking  incidence  (Emery,  Herrig,  Erwin,  and  Felix, 
1958)  are 

(O.5  +  5.O0<L,)  [3.31] 


'"stall 


=  Kef  +1'5 


10  + 


150 


and 


l  choke  l ref 


10- 


0) 


450 


(0.5  +  5.0fe),,J. 


[3.32] 


26 


The  equations  for  deviation  due  to  incidence  were  constructed  by  Hearsey  using 
existing  cascade  data  (Horlock,  1958).  These  equations  are  divided  into  four  parts, 
depending  on  the  actual  incidence  angle.  They  are 


from  Figure  3.10. 

3.3  Data  Reduction  Techniques 

Compressor  performance  characteristics  enable  a  stage-by-stage  analysis 
technique  to  describe  the  changes  in  the  flow  field  across  a  stage.  Performance 
characteristics  are  normally  specified  in  the  form  of  pressure  ratio  and  temperature 
ratio,  or  pressure  ratio  and  compressor  or  stage  efficiency,  as  a  function  of  a  mass 
flow  coefficient.  The  method  used  in  this  study  is  the  total  pressure  ratio  (XFP)  and 
total  temperature  ratio  (XPT)  as  a  function  of  corrected  mass  flow  {(/)).  Several  other 
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methods  for  calculating  total  pressure  and  total  temperature  characteristics  are 
given  in  Appendix  D. 

Because  the  data  provided  by  the  SLCC  is  a  radial  distribution  of  quantities 
defined  on  streamlines,  a  data  reduction  technique  is  required  to  generate  1-D 
performance  characteristics.  To  satisfy  the  first  law  of  thermodynamics,  a  mass- 
average  technique  must  be  used  for  total  temperature  reduction.  However,  several 
methods  are  available  for  total  pressure  ratio  averaging.  Oates  (1988)  offers  three 
different  averaging  techniques  for  total  pressure  ratio.  These  three  techniques  are 
mass-average,  stream  thrust-average,  and  continuity-average.  The  mass-average 
technique  was  chosen  for  this  study  because  of  its  simplicity.  The  mass-averaging 
technique  is  explained  next. 

The  mass  flow  between  two  streamlines  is  calculated  in  the  SLCC.  A 
linear  interpolation  between  the  characteristics  on  the  streamline  that  defines  the 
top  of  the  streamtube  and  the  streamline  that  defines  the  bottom  of  a  streamtube  is 
performed  to  find  the  characteristic  (PRt,n,  TRt,n)  in  the  center  of  the  streamtube. 
Equations  [3.37]  and  [3.38]  show  the  reduction  used  to  get  a  1-D  representation  for 
pressure  and  temperature  characteristics, 

p-37) 
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In  Equations  [3.37]  and  [3.38],  the  n  index  is  the  number  of  streamtubes  and  mt  is 
the  total  mass  flow  through  the  machine.  Figure  3.11  illustrates  how  this  reduction 
is  done. 
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4.0  CALIBRATION  AND  RESULTS  FOR  NASA 

ROTOR  IB 

To  develop  the  method  for  generating  compressor  performance 
characteristics,  a  calibration  of  the  SLCC  for  one  machine  was  required  to 
determine  the  predictive  behavior  of  the  loss  and  deviation  correlations.  The 
machine  chosen  to  do  the  initial  calibration  was  the  NASA  Rotor  IB  compressor 
that  operated  in  the  General  Electric  House  Compressor  Test  Facility,  Lynn, 
Massachusetts,  as  described  in  Section  4.1.  The  following  sections  will  discuss 
the  characteristics  of  NASA  Rotor  IB  and  its  test  setup,  the  calibration  technique 
of  the  SLCC,  and  the  results  of  the  calibration  for  Rotor  IB. 

4.1  Experimental  Data  for  NASA  Rotor  IB 

NASA  Rotor  IB  (Seyler  and  Gostolow,  1967)  was  chosen  as  the  initial 
calibration  machine  for  this  approach  for  several  reasons.  Rotor  IB  is  a  single 
blade  row  machine  with  no  inlet  or  exit  guide  vanes  to  alter  the  flow.  The  data 
provided  in  the  report  consists  of  extensive  radial  distributions  of  several  flow 
quantities  to  compare  with  the  SLCC  predicted  results.  With  only  one  blade  row  to 
calibrate,  the  effects  of  the  calibration  quantities  would  be  easier  to  discern. 
NASA  Rotor  IB  is  a  transonic  machine  that  has  a  combination  of  double-circular- 
arc  (DCA)  and  multiple-circular-arc  (MCA)  blade  shapes.  These  qualities  are 
similar  to  those  of  modem  military  fans. 
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4.1.1  Characteristics  of  NASA  Rotor  IB 

NASA  Rotor  IB  was  designed  in  the  mid-1960’s  to  advance  the  level  of 
knowledge  of  high  speed,  high  stage  loading  turbomachines.  With  the 
development  of  a  workable  titanium  alloy,  blades  were  strong  enough  to  have  tip 
speeds  beyond  1400  feet  per  second,  where  shock  losses,  due  to  high  relative 
Mach  numbers,  are  important. 

The  blades  of  Rotor  IB  were  designed  with  a  DC  A  profile  from  the  hub  to 
approximately  60-percent  span.  The  DCA  profile  is  composed  of  circular-arc 
upper  and  lower  surfaces.  The  arc  for  each  surface  is  drawn  between  the  point  of 
maximum  thickness  at  mid-chord  and  the  tangent  to  the  circles  of  the  leading-  and 
trailing-edge  radii  (Lieblein,  1965).  A  MCA  camber  line  was  used  from 
60-percent  span  to  the  tip.  The  MCA  profile  is  similar  to  a  DCA  profile  except 
that  the  two  circular  arcs  are  mutually  tangent  at  the  point  directly  across  the  flow 
passage  from  the  leading  edge  of  the  adjacent  blade  that  forms  the  other  side  of  the 
flow  passage,  instead  of  at  the  mid-chord  (Seyler  and  Gostolow,  1967).  The  rotor 
tip  diameter  was  36.5  inches  with  a  hub-tip  ratio  of  0.50,  a  tip  solidity  of  1.30,  and 
a  hub  solidity  of  about  2.45.  A  mid-span  damper  was  employed  to  maintain  the 
structural  integrity  of  each  blade.  The  rotor  had  a  medium  aspect-ratio,  a  tip  speed 
of  1400  feet  per  second,  and  a  diffusion  factor  of  0.35  at  the  tip.  The  machine  had 
a  total  of  44  blades  in  the  single  blade  row. 
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The  design  point  for  Rotor  IB  was  a  flow  of  215.49  pounds  per  second  at 
a  total  pressure  ratio  of  1.60  and  adiabatic  efficiency  of  0.858.  The  machine 
actually  achieved  a  total  pressure  ratio  of  1.638  and  adiabatic  efficiency  of  0.895  at 
a  flow  of  219.2  pounds  per  second. 

4.1.2  Test  Facility 

Rotor  IB  was  tested  at  General  Electric’s  House  Compressor  Test  Facility 
in  Lynn,  Massachusetts.  The  test  rotor  drew  air  from  the  atmosphere  through  two 
banks  of  filters.  The  first  filter  bank  was  designed  to  remove  22  percent  of  the 
particles  in  the  air  larger  than  abut  three  microns.  The  second  filter  was  designed 
to  remove  90-  to  95-  percent  of  the  remaining  particles  of  the  same  size.  The  air 
then  went  through  a  coarse  wire  foreign  object  damage  (FOD)  screen,  a  bell- 
mouth,  and  a  flow  straightener.  An  area  contraction  of  2.24  occurred  between  the 
flow  straightener  exit  and  the  inlet  face  of  the  test  rotor.  Outlet  guide  vanes 
(OGV)  far  downstream  of  the  rotor  were  used  to  de-swirl  the  flow  leaving  the  test 
facility.  The  exit  flow  was  then  split  into  two  flows,  each  of  which  was  passed 
through  a  venturi  for  flow  measurement,  and  then  exhausted  to  atmosphere.  The 
test  rotor  was  powered  by  a  high-pressure,  non-condensing  steam  turbine  rated  at 
15,000  horsepower. 

4.1.3  Instrumentation 

Inlet  conditions  to  Rotor  IB  were  measured  using  24  thermocouples 
distributed  on  the  inlet  screen  and  6  pitot-static  rakes  of  seven  elements  each 
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placed  at  centers  of  equal  annulus  areas  located  14  inches  downstream  of  the  flow 
straightener.  Blade  element  data  was  recorded  at  five  radial  positions  using 
transverse  probes.  These  measurements  were  taken  at  10-,  30-,  50-,  70-,  and  90- 
percent  of  the  annulus  height.  Immersion  of  the  transverse  probes  at  other 
instrumentation  points  were  established  to  correspond  to  the  radial  location  at 
which  the  design  streamlines  would  pass  the  instrumentation  planes.  Several  static 
pressure  taps  were  located  on  the  casing  and  hub  throughout  the  flow  path.  Static 
pressure  was  measured  at  only  one  circumferential  location,  except  at  the  axial 
stations  where  the  traverse  probes  or  rakes  were  located.  There,  the  static  pressure 
was  measured  at  more  than  one  circumferential  location.  The  thermocouple  rakes 
and  static  pressure  wedges  were  calibrated  for  Mach  number  effects  and  were 
sufficiently  insensitive  to  small  pitch  and  yaw  effects.  The  test  configuration  and 
instrument  location  can  be  seen  in  Figure  4.1. 

4.2  Calibration  Technique  for  NASA  Rotor  IB 

Four  calibration  quantities  were  chosen  by  the  author  as  the  means  to 
calibrate  the  SLCC  to  Rotor  IB  performance.  These  were  radial  variations  of 
additional  loss,  additional  deviation,  inlet  blade  blockage,  and  exit  blade  blockage. 
Additional  loss  and  deviation  parameters  correct  the  predicted  loss  and  deviation 
to  experimental  data.  The  Tioss  term  was  used  as  suggested  in  Section  3.2.1.  The 
Hioss  term  mentioned  in  Section  3.2.1  wasn’t  used  because  it  could  not  be  adjusted 
to  correctly  match  experimental  data  at  the  hub.  With  the  predicted  loss  and 
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deviation  matching  experimental  data,  the  other  predicted  flow  field  quantities 
were  much  closer  to  experimental  data.  The  radial  distributions  of  inlet  and  exit 
blade  blockages  were  used  to  achieve  the  best  match  between  predicted  blade  row 
effects  on  the  flow  and  the  experimental  data. 

Because  no  experimental  error  was  available  in  the  experimental  data  for 
Rotor  IB,  it  was  speculated  that  if  the  SLCC  results  were  within  one-percent  of 
experimental  data,  it  would  be  within  the  experimental  data  error  band.  Therefore, 
the  goal  in  the  calibration  effort  was  to  determine  the  exit  total  pressure  and  exit 
total  temperature  to  within  one-percent  of  the  data. 

The  results  of  each  calibration  step  are  presented  in  this  section  to  show 
how  each  step  in  the  calibration  effort  was  an  improvement  over  the  previous  step. 
The  quantities  discussed  are  loss,  deviation,  exit  absolute  Mach  number,  exit 
relative  Mach  number,  exit  absolute  swirl  velocity,  exit  relative  swirl  velocity,  exit 
total  pressure,  exit  static  pressure,  exit  total  temperature,  and  exit  static 
temperature.  These  parameters  were  chosen  because  they  give  a  complete 
description  of  the  change  in  the  flow  field  through  a  bladed  region  in  terms  of 
changes  occurring  in  the  absolute,  relative,  total,  and  static  reference  frames. 
Therefore,  they  give  a  good  overall  picture  of  machine  performance. 

The  three  levels  of  calibration  in  the  calibration  effort  are  denoted 
uncalibrated,  partially  calibrated,  and  calibrated.  The  results  of  the  SLCC  with  no 
calibration  are  presented  as  a  baseline  (uncalibrated)  to  see  how  well  the 
uncorrected  correlations  and  the  SLCC  performed  to  predict  blade  row 
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performance.  Figure  4.2(a)  illustrates  how  this  method  of  calibration  was 
conducted.  The  only  input  was  geometry.  The  SLCC  was  then  allowed  to  run 
with  only  the  correlations  used  to  predict  the  results. 

Next,  because  experimental  loss  and  deviation  were  available  from  the 
Rotor  IB  data,  the  SLCC  was  altered  to  accept  loss  and  deviation  as  user  inputs. 
The  correlations  were  still  allowed  to  calculate  a  loss  and  deviation  at  each 
streamline.  The  difference  between  the  loss  and  deviation  predicted  by  the 
correlations  and  experimental  loss  and  deviation  was  included  as  add-loss  and  add- 
deviation.  Performance  results  with  only  the  loss  and  deviation  corrections 
utilized  are  referred  to  as  the  partially  calibrated  results.  Figure  4.2(b)  illustrates 
how  this  calibration  was  conducted.  Geometry,  experimental  loss,  and 
experimental  deviation  were  input  to  the  SLCC.  The  results  of  this  run  yielded 
add-loss  and  add-deviation  as  described  above.  The  geometry,  add-loss,  and  add- 
deviation  were  then  input  into  the  SLCC  to  yield  the  final  results. 

The  last  calibration  step  was  the  fully  calibrated  step.  For  this  step, 
geometry,  add-loss,  add-deviation,  inlet  blade  blockage  and  exit  blade  blockage 
were  input  to  the  SLCC.  The  SLCC  results  were  then  compared  to  experimental 
results.  If  the  results  were  within  the  predetermined  one-percent  tolerance,  the 
results  were  denoted  as  final  results.  If  the  results  were  outside  of  the  one-percent 
tolerance,  the  blockage  terms  were  adjusted  and  the  SLCC  was  initiated  again. 
This  was  an  iterative  process  to  adjust  the  blockage  terms.  The  process  for  the 
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fully  calibrated  step  is  illustrated  in  Figure  4.2(c).  The  performance  results  with 
all  four  calibration  parameters  included  are  referred  to  as  the  calibrated  results. 
4.2.1  Description  of  Calibration  Points  and  Figures 

The  results  of  the  calibration  are  shown  in  Figures  4.3  through  4.12.  Each 
of  the  three  calibration  steps  is  shown  in  comparison  to  experimental  data.  The 
experimental  data  is  displayed  with  black  diamonds,  the  uncalibrated  results  are 
represented  with  a  green  dot-dashed  line,  the  partially  calibrated  results  are  shown 
with  a  dashed  red  line,  and  the  fully  calibrated  results  are  shown  with  a  solid  blue 
line.  Table  4.1  summarizes  the  percent  difference  between  experimental  data  and 
the  results  from  each  of  the  three  calibration  steps.  The  loss  results  in  Table  4.1 
are  given  in  absolute  error  instead  of  percent  error  because  the  values  are  so  close 
to  zero  that  percent  differences  would  give  an  unclear  representation  of  the  results. 

Although  the  results  were  predicted  for  50-,  70-,  90-,  and  100-percent 
design  speeds  and  numerous  mass  flows  along  each  speed  line,  only  one  speed  and 
flow  point  is  shown  for  comparison.  The  point  chosen  is  on  the  100-percent  speed 
line  with  a  corrected  mass  flow  of  216.  31  lbm/s.  This  point  was  chosen  because 
the  100-percent  speed  line  represented  the  most  difficult  speed  to  calibrate 
(apparently  due  to  the  transonic  flows  relative  to  the  blade  causing  increase  shock 
losses).  The  chosen  point  is  also  in  the  middle  of  the  speed  line,  and  considered  to 
be  typical  of  most  points  along  the  100-percent  speed  line.  Predictions  at  other 
speeds  were  typically  better  than  the  ones  presented  in  this  section. 
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4.2.2  Discussion  of  Calibration  Results 

The  uncalibrated  results  for  the  relative  total  pressure  loss  coefficient, 
shown  in  Figure  4.3,  show  that  the  correlations  significantly  over-predict  the  loss 
around  the  hub.  Rotor  IB  has  a  hub  solidity  approaching  2.5.  The  correlations  are 
only  reliable  up  to  a  solidity  of  about  2.0.  Beyond  2.0,  the  correlations  become 
unstable.  This  is  why  the  loss  was  over-predicted  at  the  hub.  However,  the 
uncalibrated  predictions  are  closer  toward  the  tip.  For  both  the  partial  calibration 
and  full  calibration  steps,  the  loss  is  specified.  The  loss  is  still  in  error  by  a  small 
amount  (0.013  and  0.014,  respectively)  as  opposed  to  -0.15  for  the  uncalibrated 
prediction.  This  error  in  the  partially  and  full  calibration  prediction  is  partly  due  to 
interpolation  and  extrapolation  of  experimental  data.  Interpolation  and 
extrapolation  are  required  since  experimental  loss  data  is  given  at  only  five  radial 
locations  across  the  flow.  Losses  at  locations  of  the  streamlines  predicted  by  the 
SLCC  were  interpolated  or  extrapolated  from  the  five  experimental  points.  The 
partially  calibrated  and  calibrated  results  in  Figure  4.3  lie  nearly  on  top  of  each 
other. 

Figure  4.4  shows  that  the  correlations  under-predict  the  deviation  at  the 
hub,  and  again,  do  better  near  the  tip  for  the  uncalibrated  step.  These  results  are 
somewhat  expected  because  the  correlations,  for  the  most  part,  were  developed 
from  linear  cascade  theory  and  experience  and  the  flow  phenomena  occurring  at 
hub  and  tip  regions  are  highly  three-dimensional.  For  the  partially  calibrated  and 
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fully  calibrated  results,  where  the  deviation  is  specified,  the  predictions  shows  a 
30-percent  improvement  over  the  uncalibrated  results.  The  predicted  results  differ, 
in  part,  from  the  experimental  data  also  because  of  required  interpolation  or 
extrapolation  of  theoretical  results  from  experimental  data. 

Exit  absolute  Mach  number,  shown  in  Figure  4.5,  is  fairly  constant  across 
the  radius  of  the  machine.  The  uncalibrated  values  differ  from  experiment  by  a 
maximum  of  -3.87-percent  in  the  hub  region.  With  the  partial  calibration,  an 
improvement  of  1.73-percent  is  gained  in  the  prediction  over  the  uncalibrated  run. 
The  complete  calibration  prediction  gives  a  slightly  better  improvement  to  a 
-1.93-percent  difference  from  experimental  data. 

The  exit  relative  Mach  numbers,  shown  in  Figure  4.6,  differ  from 
experiment  by  9.88-percent  at  the  compressor  hub  for  the  uncalibrated  case.  In  the 
tip  region,  the  values  more  closely  follow  the  actual  data.  The  partially  calibrated 
prediction  improves  by  more  than  50-percent  over  the  uncalibrated  error.  The 
fully  calibrated  predictions  exhibit  an  error  of  3.79-percent  for  the  exit  relative 
Mach  number. 

The  exit  absolute  swirl  velocity  is  one  of  the  quantities  needed  by  the 
SLCC  to  step  across  the  bladed  region.  The  radial  distribution  for  it  is  shown  in 
Figure  4.7.  The  prediction  of  the  exit  absolute  swirl  velocity  for  the  uncalibrated 
case  differs  a  substantial  -11.05-percent  from  experiment  in  the  hub  region.  This 
quantity  is  directly  dependant  on  predicted  deviation,  which  differs  significantly 
from  experiment  at  the  hub.  With  loss  and  deviation  specified,  the  error  for  the 
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partially  calibrated  case  decreases  to  8.06-percent  at  the  tip.  With  the  radial  inlet 
and  exit  blockage  included  in  the  full  calibration,  the  error  comes  down  to  a  more 
reasonable  -5.32-percent. 

The  exit  relative  swirl  velocity  is  shown  in  Figure  4.8.  For  the  uncalibrated 
case,  the  predicted  values  differ  from  experiment  by  a  large  amount,  22.70-percent 
at  the  hub.  The  partial  and  calibrated  results  are  fairly  close,  with  the  calibrated 
results  reducing  the  percent  error  by  almost  a  factor  of  seven  to  3.25-percent. 

The  exit  absolute  total  pressure  is  shown  in  Figure  4.9.  These  values  are 
one  of  the  key  factors  in  determining  performance  characteristics  and  are  another 
of  the  closure  quantities.  Because  loss  is  used  directly  in  calculating  this 
parameter,  the  exit  total  pressure  is  expected  to  differ  from  experimental  data.  The 
error  in  the  uncalibrated  run  is  2. 81 -percent.  However,  with  loss  and  deviation 
corrected,  the  error  in  the  predicted  exit  total  pressure  is  worse,  with  ati  error  of 
—3.08-percent  when  compared  to  experiment.  The  results  match  much  better  at  the 
hub,  but  miss  more  at  the  tip.  The  author  interprets  these  results  to  mean  that  the 
inlet  and  exit  radial  blockages  are  additional  important  factors  in  determining  exit 
total  pressure.  With  a  complete  calibration,  the  percent  error  for  exit  total  pressure 
compared  to  experimental  data  is  -0.89-percent.  This  is  within  the  1.0-percent 
goal  set  earlier. 

Figure  4.10  shows  the  exit  static  pressure.  The  uncalibrated  run  differs  a 
maximum  of  3.35-percent  at  the  hub.  The  partial  and  calibrated  runs  lie  nearly  on 
top  of  each  other  until  the  radius  approaches  16-inches.  The  partial  results  then 
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over-predict  the  data  by  over  2.0-percent.  The  complete  calibration  brings  the 
maximum  error  to  -0.56-percent  at  the  radial  location  of  about  16-inches. 

The  last  of  the  closure  quantities  is  the  exit  absolute  total  temperature, 
shown  in  Figure  4.1 1.  The  uncalibrated  SLCC  predictions  differ  by  a  maximum  of 
-1.36-percent.  The  predicted  temperature  is  high  at  the  hub,  low  in  the  middle, 
and  high  at  the  tip.  The  partial  calibration  reduces  the  error  to  less  than 
1.0-percent.  The  final  calibrated  results  for  the  exit  total  temperature  differ  from 
the  experimental  data  by  a  maximum  of  0.18-percent. 

Exit  static  temperature  is  plotted  in  Figure  4.12.  The  uncalibrated  results 
followed  the  same  trends  as  the  uncalibrated  total  temperature  results.  The 
maximum  error  for  the  exit  static  temperature  for  the  uncalibrated  case  were 
within  the  one-percent  goal.  The  partial  and  calibrated  results  showed  slight 
improvements  with  a  final  percent  error  of  only  0.10-percent. 

In  summary,  only  the  value  of  the  exit  static  temperature  came  in  within 
the  desired  one-percent  accuracy  for  the  uncalibrated  parameters.  The  results  of 
the  partial  calibration  typically  showed  significant  improvements  over  the  baseline 
correlation  predictions  alone.  The  final  calibration  showed  small  improvements 
over  the  partial  calibration.  Not  all  values  tracked  were  less  than  the  one-percent 
goal.  However,  the  exit  total  pressure  and  exit  total  temperature  are  the  main 
parameters  used  in  compressor  performance  characteristic  prediction  and  they 
were  both  within  the  one-percent  goal.  Again,  the  calibration  point  discussed  in 
this  section  was  for  a  high  speed  point  where  the  correlations  are  expected  to  have 
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the  most  problems  (because  of  transonic  flow  conditions).  These  results  are 
typical  for  the  100-percent  speed  line.  At  lower  speeds,  the  error  between 
experimental  and  calibrated  data  points  were  typically  lower  than  the  point 
presented  here. 

4.3  Calibration  Trends 

The  above  calibration  process  was  described  for  one  speed  and  flow  point 
only.  However,  this  calibration  method  was  also  carried  out  for  50-,  70-,  90-,  and 
100-percent  speeds  at  several  mass  flows  from  the  stall  line  to  the  choking  region 
for  each  speed.  Several  trends  became  apparent  throughout  this  process.  These 
trends  are  discussed  in  this  section. 

4.3.1  Trends  Along  the  100-percent  Speed  Line 

A  radial  distribution  of  additional  loss  (add-loss)  was  included  to  calibrate 
the  loss  correlations  to  match  experimental  data.  Add-loss  is  loss  added  (or 
subtracted)  to  the  loss  calculated  by  the  correlations.  The  calculation  of  the  add- 
loss  term  is  described  in  Section  4.2.  Add-loss  was  included  at  every  experimental 
point  on  the  map.  The  radial  distribution  of  add-loss  for  the  point  at  100-percent 
speed  is  shown  in  Figure  4.13.  There  are  seven  different  mass  flows  on  this  figure 
ranging  from  204.79  lbm/s  to  221.81  lbm/s.  An  eighth  curve  is  also  shown  on  the 
figure.  It  is  a  solid  heavy  black  line  that  represents  a  least-squares  curve  fit  to  the 
data.  As  shown  in  Figure  4.3,  the  loss  at  the  hub  is  over-predicted  with  the 
baseline  correlations.  Figure  4.13  clearly  shows  how  the  loss  is  subtracted  at  the 
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hub  to  account  for  this  effect.  From  the  mid-span  to  the  tip,  the  additional  loss 
term  remains  fairly  constant.  The  higher  mass  flows  do  cause  the  additional  loss 
to  decrease  toward  the  tip.  The  overall  fit  curve  has  a  parabolic  shape  with 
negative  values  at  the  hub,  increasing  to  positive  values  past  the  mid-span,  and 
then  decrease  into  the  negative  region  again  at  the  tip. 

Additional  Deviation  (add-deviation)  was  the  second  quantity  to  be 
calibrated.  Add-deviation  is  deviation  added  (or  subtracted)  to  the  deviation 
calculated  by  the  correlations.  The  calculation  of  the  add-deviation  term  is 
described  in  Section  4.2.  Figure  4.14  shows  the  add-deviation  data  for  100-percent 
speed.  The  curves  show  lines  of  add-deviation  for  the  same  range  of  the  mass 
flows  investigated  for  the  add-loss  factor.  Like  the  loss  correlations,  the  deviation 
correlations  apparently  have  the  most  trouble  at  the  hub.  The  add-deviation  curves 
also  have  a  parabolic  shape,  with  the  largest  values  positive  at  the  hub,  decreasing 
to  almost  zero  at  the  mid-span,  and  increasing  again  slightly  at  the  tip.  The  lower 
the  mass  flow  the  less  additional  deviation  needs  to  be  added. 

Figures  4.15  and  4.16  are  plots  of  the  radial  inlet  and  exit  blade  blockage 
terms.  These  terms  were  added  to  adjust  the  physical  area  blockage  of  the  blades 
in  the  flow  field.  For  calibration  purposes,  the  values  were  allowed  to  be  either 
positive  or  negative.  The  positive  values  are  easily  explained  as  an  actual  physical 
blockage  of  the  flow  caused  by  a  blade  and  its  boundary-layers.  The  negative 
values,  however,  are  not  inherently  obvious.  The  correlations  used  in  the  SLCC 
use  the  solidity  of  the  bladed  regions.  This  solidity  increases  toward  the  hub  of  the 
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compressor.  The  correlations  used  were  only  reliable  to  a  solidity  of  about  two. 
Near  the  hub,  the  value  of  solidity  for  Rotor  IB  exceeds  the  value  of  2.0  to  a  value 
of  nearly  2.5.  This  causes  the  correlations  predict  unreasonable  amount  of  loss  at 
the  hub.  The  negative  values  of  blockage  help  to  adjust  for  the  effects  of  the  high 
values  of  solidity  and  other  3-D  flow  effects  that  cause  problems  in  the 
correlations.  Static  pressure  and  numerous  velocity  values  were  identified  as  the 
quantities  most  affected  by  blockage.  An  increase  in  blockage  caused  an  increase 
in  the  local  velocity  and  a  decrease  in  the  local  static  pressure.  The  inlet  blockage 
of  the  blade  row  had  little  effect  on  the  overall  flow  field.  However,  the  entire 
flow  field  was  very  dependant  on  the  exit  blockage. 

Figure  4.15  shows  the  inlet  blockage  used  to  completely  calibrate  the 
SLCC  for  Rotor  IB.  Inlet  blockage  affected  the  inlet  static  pressure  and 
temperature  and  velocity  quantities  but  had  little  effect  on  the  exit  conditions.  The 
inlet  blockage  at  the  hub  is  negative  and  increases  to  positive  just  before  the  mid¬ 
span.  From  the  mid-span  to  the  tip,  the  value  of  inlet  blockage  remains  fairly 
constant  at  about  4.0-percent. 

The  exit  blockage  (Figure  4.16)  was  used  to  adjust  the  exit  static  pressure 
and  the  velocities.  This  calibration  quantity  seemed  to  be  more  important  to 
overall  calibration  than  the  inlet  blockage.  The  variation  of  exit  blockage  affected 
both  inlet  and  exit  conditions  significantly.  The  shape  of  this  curve  is  very  similar 
to  the  inlet  blockage  except  it  increases  more  significantly  at  the  tip.  It  is  negative 
at  the  hub  and  increases  to  about  zero,  and  remains  around  zero  from  about 
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20-percent  to  80-percent  span.  Beyond  80-percent  span,  the  blockage  increases 
again.  This  increase  in  blockage  decreases  the  predicted  static  pressure  to  match 
the  experimentally  determined  static  pressure.  The  fit  line  is  again  a  least-square 
curve  fit  to  all  of  the  different  mass  flows  along  the  given  speed  line. 

4.3.2  Trends  Over  the  Entire  Machine 

The  calibration  trends  presented  in  the  preceding  section  were  duplicated 
for  the  entire  machine.  The  trends  presented  in  the  previous  section  were  fairly 
typical  of  the  results  seen  at  other  speeds.  The  calibration  results  presented  in 
Figures  4.17  through  4.20  are  the  fit  lines  for  each  speed  and  calibration  quantity. 
To  make  the  results  easily  applicable  to  another  machine,  the  results  are  presented 
as  a  function  of  percent  span  instead  of  radius. 

Figure  4.17  shows  the  radial  distribution  of  the  averaged  add-loss  for 
different  speeds.  Each  curve  has  a  similar  parabolic  shape,  with  the  values  largely 
negative  around  the  hub,  increasing  to  the  mid-span,  and  then  decreasing  to  the  tip. 
The  magnitude  of  add-loss  decreases  with  decreasing  machine  speed.  This  is 
expected  because  most  of  the  correlations  used  were  originally  developed  for  low 
speed  machines.  Around  the  mid-span,  the  curves  for  all  the  speeds  converge 
around  the  same  values. 

The  additional  deviation  is  presented  in  Figure  4.18.  It  too  has  a  parabolic 
shape,  decreasing  toward  the  mid-span  and  then  increasing  to  the  tip.  The  curves 
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again  converge  to  similar  values  at  the  mid-span.  The  additional  deviation  needed 
typically  decreases  with  decreasing  speed. 

Radial  inlet  blockage  as  a  function  of  percent  span  is  shown  in  Figure  4.19. 
All  of  the  curves  have  similar  shapes  and  similar  magnitudes.  They  are  largely 
negative  at  the  hub,  increasing  toward  the  mid-span  (where  they  are  fairly 
constant),  and  increasing  slightly  at  the  tip. 

The  final  calibration  curves,  exit  blockage,  are  seen  in  Figure  4.20.  Like 
the  inlet  blockage,  the  exit  blockage  curves  are  very  similar  in  shape  and 
magnitude.  Each  curve  tends  to  have  a  negative  value  at  the  hub  and  increase  to 
almost  zero  at  the  mid-  span.  From  mid-span  to  tip,  they  increase  at  a  much  more 
aggressive  rate  than  the  inlet  blockage.  Each  of  the  curves  tends  to  lay  almost  on 
top  of  the  others  throughout  the  span  except  for  the  70-percent  curve. 

4.4  Compressor  Maps 

With  the  predicted  radial  flow  field  of  Rotor  IB  calibrated,  the  computed 
flow  fields  were  processed  through  a  mass-averaged  data  reduction  routine  to 
reduce  it  to  a  1-D  representation.  With  the  data  in  a  1-D  form,  it  was  possible  to 
generate  1-D  compressor  characteristics.  The  characteristics  presented  are  the 
total  pressure  ratio  (PR)  and  total  temperature  ratio  (TR)  versus  corrected  mass 
flow  for  different  speeds.  These  results  are  shown  in  Figures  4.21  and  4.22, 
respectively. 
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The  pressure  characteristic  shown  for  50-,  70-,  90-,  and  100-percent  design 
corrected  speed  are  plotted  in  Figure  4.21.  Corrected  mass  flow  (<p)  is  on  the 
abscissa  and  PR  (xPp)  on  the  ordinate.  The  calibrated  data  matches  the 
experimental  data  closely.  The  maximum  percent  error  for  every  point  on  the  50-, 
70-,  and  90-percent  speed  lines  is  below  one-percent.  On  the  highest  speed  line, 
100-percent,  all  points  are  below  one-percent  except  for  the  highest  mass  flow 
point,  which  is  in  the  choking  region.  The  maximum  percent  error  at  this  point  is 
2.63-percent.  This  point  could  not  be  completely  calibrated.  Any  further 
adjustment  to  the  calibration  quantities  would  not  allow  the  SLCC  to  converge  to  a 
solution.  In  general,  the  highest  percent  difference  on  a  speed  line  is  the  point 
closest  to  choking.  This  is  to  be  expected  since  the  streamline  curvature  method 
has  difficulty  with  choking  (Cumpsty,  1989).  The  error  also  typically  increases  as 
the  speed  increases. 

Figure  4.22  shows  the  temperature  characteristics  for  50-,  70-,  90-,  and 
100-percent  corrected  speed.  The  abscissa  is  again  (f)  but  the  ordinate  is  TR  (xFt). 
These  characteristics  have  the  same  general  shape  as  the  pressure  characteristics, 
however  the  magnitudes  are  not  quite  as  large.  The  errors  followed  trends  similar 
to  those  for  the  pressure  characteristic:  errors  tend  to  increase  toward  the  choking 
region  and  as  speed  increases.  The  percent  error  at  each  point  is  under 
one-percent.  The  largest  error,  0.67-percent  is  again  at  the  choking  point  on  the 
100-percent  speed  line. 
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The  next  two  figures  show  the  2-D  characteristics.  The  characteristics  are 
in  the  same  form  as  the  ones  presented  in  Figures  4.21  and  4.22.  However,  the 
results  are  no  longer  a  1-D  representation,  but  are  broken  into  hub,  mid-line,  and 
tip  regions.  The  percent  span  from  the  hub  are  10-,  50-,  and  90-percent.  On  the 
figures,  diamonds  represent  the  experimental  data,  a  red  dot-dash  line  represents 
the  hub  results,  a  solid  blue  line  represents  the  mid-line,  and  a  green  dotted  line 
represents  the  tip  results. 

Figure  4.23  shows  the  pressure  characteristics.  At  lower  speeds,  50-  and 
70-percent,  the  characteristic  curves  lie  close  to  each  other  and  are  relatively  flat. 
As  the  mass  flow  increases  at  the  lower  speeds,  they  diverge  some  but  are  still 
relatively  close  and  flat.  This  means  that  at  the  lower  speeds,  there  is  little 
variation  in  work  done  across  the  span  of  the  blade.  At  the  higher  speeds,  the  90- 
and  100-percent  lines,  the  character  of  the  curves  changes.  The  tip  region  begins 
to  do  more  work  at  lower  mass  flows.  As  the  mass  flow  approaches  choked  flow 
conditions,  the  characteristics  converge.  The  characteristic  curves  at  the  hub  tend 
to  be  flatter,  while  the  tip  drops  off  faster  as  the  flow  moves  toward  choke. 

The  percent  error  of  pressure  characteristics  follows  the  same  trends  as  for 
the  1-D  representation.  The  maximum  percent  difference  occurs  at  the  points 
closest  to  choke  for  100-percent  speed  at  the  tip.  From  looking  at  the  2-D  pressure 
characteristics,  the  tip  area  seems  to  be  more  into  the  choking  region  than  the  mid¬ 
line  or  hub.  Since  the  streamline  curvature  method  has  difficulty  with  choking, 
this  point  can  be  expected  to  be  less  accurate.  The  percent  errors  in  all  of  the 
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pressure  characteristics  at  the  hub  are  less  than  0.9-percent.  The  maximum  percent 
error  for  the  mid-line  is  2.35-percent  and  5.06-percent  for  the  tip. 

The  temperature  characteristics,  shown  in  Figure  4.24,  also  follow  the  same 
trends  as  the  pressure  characteristics.  At  lower  speeds,  the  curves  are  flatter  and 
lie  closer  together.  As  speed  increases,  the  hub,  mid-line,  and  tip  take  on  a 
different  character.  They  again  converge  as  the  flow  tends  toward  choked 
conditions.  The  maximum  error  for  the  entire  temperature  characteristics  map  is 
again  at  the  100-percent  speed  line  at  the  point  closest  to  choke  in  the  tip  region. 
The  error  at  this  point  is  1.37-percent.  The  error  for  all  other  points  on  the  map  is 
less  than  0.60-percent. 

These  results  show  that  with  the  relative  total  pressure  coefficient, 
deviation  angle,  and  inlet  and  exit  blockage  calibrated  in  the  SLCC,  the  flow  field 
is  approximated  closely.  With  a  calibrated  flow  field,  the  1-D  and  2-D 
characteristics  are  shown  to  be  in  good  agreement  with  the  experimental  data. 
Therefore,  a  technique  to  predict  stage-by-stage,  pre-stall  compressor  performance 
characteristics  has  been  developed. 
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5.0  PREDICTION  OF  NASA  STAGE  35 
COMPRESSOR  MAPS 

With  NASA  Rotor  IB  calibrated,  it  was  desirable  to  use  the  knowledge 
gained  from  that  calibration  to  predict  the  performance  of  some  other  machine. 
An  investigation  was  initiated  to  determine  if  the  calibration  trends  discovered  for 
NASA  Rotor  IB  are  useful  for  calibrating  a  very  different  machine  such  as  NASA 
Stage  35.  The  following  section  describes  NASA  Stage  35  and  shows  the  radial 
and  overall  predictions  and  comparisons  of  the  pre-stall  compressor  performance 
maps. 

5.1  Experimental  Data  for  NASA  Stage  35 

NASA  Stage  35  (Reid  and  Moore,  1978)  was  chosen  as  the  machine  to 
predict  performance  by  applying  calibration  trends  gained  by  analysis  of  NASA 
Rotor  IB.  Stage  35  is  a  single  stage,  two  blade  row  machine  designed  in  the  mid- 
1970’s.  It  was  chosen  because  it  has  highly  loaded,  low  aspect  ratio  blades  that  are 
typical  in  modern  military  high  pressure  compressor  cores.  There  is  also 
experimental  data  available  to  compare  with  the  prediction.  Rotor  IB  is  designed 
to  be  one  of  the  rotors  of  a  low  speed  compressor  fan  stage,  while  Stage  35  is 
designed  to  be  the  first  stage  of  a  core  compressor.  The  rotor  tip  speed  is  similar 
to  that  of  Rotor  IB.  However,  Rotor  IB  is  a  moderate-aspect  ratio  machine 
whereas  Stage  35  is  a  low  aspect-ratio  machine.  These  different  characteristics 
will  verify  if  similar  calibration  quantities  will  hold  for  dissimilar  machines. 
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5.1.1  Characteristics  of  NASA  Stage  35 

NASA  Stage  35  was  designed  in  the  mid- 1 970’ s  for  an  investigation  into 
the  effects  of  blade  aspect  ratio  and  stage  loading.  Four  machines  were  designed 
with  identical  rotor  and  stator  solidity,  rotational  speeds,  mass  flows,  and  flow  path 
geometry.  Only  the  rotor  aspect  ratio  and  stage  loading  were  varied.  Stage  35  has 
a  low-aspect  ratio  rotor  and  low  pressure  ratio  design. 

Both  rotor  and  stator  blades  were  designed  with  multiple-circular  arc 
(MCA)  profiles  and  low  aspect  ratios.  The  rotor  tip  diameter  is  19.9  inches  with  a 
hub-tip  ratio  of  0.7.  Because  of  the  higher  hub-tip  ratio,  no  mid-span  damper  was 
used.  The  machine  has  36  rotor  blades  and  46  stator  blades.  The  tip  solidity  of  the 
rotor  and  stator  regions  are  both  1.3.  The  hub  solidity  for  the  rotor  and  stator  is 
1.77  and  1.48,  respectively. 

The  design  point  for  Stage  35  is  a  mass  flow  of  44.53  pounds  per  second  at 
a  stage  total  pressure  ratio  of  1.82.  The  peak  rotor  and  stage  efficiency  are  0.872 
and  0.845,  respectively.  The  rotor  tip  speed  is  1493  feet  per  second  and  has  a 
maximum  rotor  diffusion  factor  of  0.48.  The  stator  diffusion  factor  is  0.34. 

5.1.2  Test  Facility 

Stage  35  was  tested  at  NASA  Lewis  Research  Center’s  Single-Stage 
Compressor  Facility  in  Cleveland,  Ohio.  A  schematic  of  the  facility  is  show  in 
Figure  5.1.  Atmospheric  air  enters  the  facility  from  the  roof  and  passes  through  a 
flow  measuring  orifice  and  into  a  plenum  chamber  upstream  of  the  test  section. 
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Air  then  passes  through  the  test  section  and  into  the  collector  and  exits  the  facility 
through  a  vacuum  exhaust  system. 

5.1.3  Instrumentation 

The  radial  distribution  of  flow  conditions  upstream  of  the  rotor  were 
determined  using  two  combination  probes  and  two  18°  wedge  probes.  The 
combination  probes  measured  total  temperature,  total  pressure,  and  flow  angle. 
The  wedge  probes  measured  static  pressure  and  flow  angle.  The  probes  were 
automatically  aligned  in  the  flow  direction  with  a  null-balancing  control  system. 
Inner  and  outer  wall  static  pressure  taps  were  located  at  the  same  axial  position  as 
the  survey  probes.  Chromel-constantan  thermocouples  were  used  to  measure 
temperature.  No  data  were  taken  between  the  rotor  and  stator  because  of  the  close 
blade  spacing.  A  radial  distribution  of  flow  conditions  at  the  stator  exit  was 
determined  using  the  two  combination  and  two  wedge  probes.  The  probes  were 
also  traversed  circumferentially.  The  layout  for  locations  of  the  combination 
probe,  wedge  probe,  and  wall  static  taps  is  shown  in  Figure  5.2.  Station  1  is  the 
station  immediately  in  front  of  the  rotor.  Station  3  is  the  station  directly  behind  the 
stator.  The  view  of  both  stations  is  looking  downstream  in  the  flow. 

The  mass  flow  was  measured  by  using  a  calibrated  thin-plate  orifice.  The 
average  of  two  Chromel-constantan  thermocouples  was  used  to  obtain  the  orifice 
temperature.  Orifice  pressures  were  measured  by  calibrated  pressure  transducers. 
An  electronic  speed  counter  and  magnetic  pickup  were  used  to  measure  rotational 
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speed.  Estimated  errors  in  the  data  based  on  inherent  inaccuracies  of  the 
instrumentation  and  recording  systems  are  shown  in  Table  5.1. 

5.2  Radial  Predictions  and  Compressor  Maps 

With  the  calibration  trends  discovered  during  the  calibration  of  Rotor  IB,  it 
was  hoped  the  same  trends  for  add-loss,  add-deviation,  inlet  blockage,  and  exit 
blockage  would  hold  true  for  Stage  35.  This  was,  however,  not  the  case.  With  the 
calibration  results  of  Rotor  IB  imposed  on  Stage  35,  the  results  were 
unsatisfactory.  Figure  5.3  shows  the  1-D  mass-averaged  stage  pressure 
characteristic  for  Stage  35  with  Rotor  IB  calibration  trends  imposed.  The  results 
are  adequate  for  the  70-percent  speed  line  but  unacceptable  for  the  90-  and 
100-percent  speed  line.  The  maximum  error  in  the  pressure  characteristic  is 
-11.10-percent  on  the  100-percent  speed  line. 

Figure  5.4  shows  the  1-D  mass-averaged  stage  temperature  characteristic 
for  Stage  35  with  Rotor  IB  calibration  trends  imposed.  Like  the  pressure 
characteristic,  the  temperature  characteristic  is  acceptable  for  the  70-percent  speed 
line.  However,  the  90-  and  100-percent  speed  line  predictions  are  in  disagreement 
with  experimental  data.  The  maximum  error  in  the  temperature  characteristic  is 
-3 .41 -percent  on  the  100-percent  speed  line. 

It  was  discovered  during  the  calibration  of  Rotor  IB  that  the  percent  error 
in  the  mass-averaged  1-D  prediction  was  smaller  than  the  prediction  of  a  radial 
distribution.  Because  of  this  fact  and  the  fact  that  the  SLCC  couldn’t  converge  to 
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a  solution  at  all  flow  points  with  Rotor  IB  trends  imposed  on  Stage  35,  the 
conclusion  was  made  that  Rotor  IB  calibration  trends  are  not  acceptable  for  Stage 
35.  Therefore,  an  investigation  was  conducted  to  explain  why  the  Rotor  IB 
calibration  trends  would  not  hold  true  for  Stage  35.  The  investigation  revealed 
that  the  machines  were  too  dissimilar  in  design  for  the  same  trends  to  be  true  for 
both.  The  major  difference  discovered  was  the  solidity  of  the  blade  rows.  Rotor 
IB  had  a  relatively  high  solidity,  especially  at  the  hub  where  it  approached  2.5. 
The  correlations  used  in  the  SLCC  were  only  reliable  up  to  a  solidity  of  about  2.0. 
The  calibration  efforts  for  Rotor  IB  were  primarily  to  overcome  the  effects  of  the 
high  solidity  at  the  hub.  Negative  blockage  values  at  the  hub  were  required 
because  of  this  high  solidity.  The  hub  solidity  for  Stage  35  peaked  at  about  1.8  for 
the  rotor. 

Another  problem  encountered  in  applying  the  correlations  to  Rotor  IB  was 
the  relatively  high  3-D  effects  at  the  hub.  The  slope  of  the  hub  wall  was 
significant.  The  slope  of  the  wall  was  larger  at  the  hub  for  Rotor  IB  than  for  Stage 
35.  SLCCs  are  also  known  to  have  problems  with  lower  hub-casing  ratio  fans 
(Cumpsty,  1989)  because  of  the  strong  curvature  in  the  meridional  plane  and  the 
larger  slope  of  the  wall  at  the  hub.  Stage  35  has  a  40.0-percent  higher  hub-casing 
ratio  and  a  30.6-percent  smaller  slope  at  the  hub  wall  than  Rotor  IB. 

A  final  major  difference  between  the  two  machines  is  the  blade  profiles. 
The  majority  of  the  Rotor  IB  blade  profile  is  the  DCA  design.  Stage  35  blades 
were  designed  with  the  MCA  design  across  the  entire  span  of  the  blades.  These 
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profile,  solidity,  and  annulus  wall  shape  differences  were  enough  to  disallow  the 
same  calibration  trends  to  be  used  on  the  different  machines.  The  Rotor  IB 
correlation  corrections  moved  loss  and  deviation  in  the  wrong  direction  for  Stage 
35.  Therefore,  it  was  decided  to  predict  Stage  35  with  the  original,  uncalibrated 
correlations.  The  correlations  and  the  SLCC  were  intended  for  use  on  machines 
with  an  annulus  geometry  like  Stage  35  more  so  than  Rotor  IB.  The  results 
introduced  in  this  section  are  for  an  uncalibrated  prediction  of  Stage  35. 

The  results  of  the  uncalibrated  SLCC  for  a  single  run  are  presented  next. 
The  results  displayed  are  for  a  case  at  100-percent  design  speed  and  a  mass  flow  of 
44.69  lbm/s.  The  results  are  broken  into  those  at  the  exit  of  the  rotor  and  those  at 
the  exit  of  the  stator.  The  maximum  difference  presented  in  the  figures  is  the 
percent  difference  between  the  experimental  data  and  the  SLCC  results.  Figures 
5.5  through  5.14  are  results  at  the  exit  of  the  rotor  and  Figures  5.15  through  5.22 
are  results  at  the  exit  of  the  stator. 

Figure  5.5  displays  the  radial  distribution  of  the  relative  total  pressure  loss 
coefficient  for  Stage  35  rotor.  These  results  follow  the  trend  of  and  nearly  match 
the  experimental  data  with  no  calibration.  Figure  5.6  shows  the  comparison  for  the 
exit  flow  angle  for  Stage  35  instead  of  the  deviation.  The  exit  flow  angles  differ 
from  the  respective  deviations  by  a  constant,  the  blade  metal  angle  at  the  exit.  The 
exit  flow  angle  results  (Figure  5.6)  deviate  from  the  experimental  data,  especially 
at  the  hub.  The  maximum  difference  between  experimental  data  and  SLCC  results 
for  the  exit  flow  angle  is  26.23-percent  at  the  hub.  One  explanation  offered  for 
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this  variance  in  exit  flow  angle  is  the  blade  shape  design.  The  correlations  were 
developed,  for  die  most  part,  on  machines  with  DCA  profiles,  whereas  Stage  35 
has  MCA  profiles.  Therefore,  the  result  was  that  the  SLCC  couldn’t  match  the 
deviation.  Because  compressor  work  is  directly  related  to  the  exit  flow  deviation 
and  deviation  poorly  predicted  at  the  hub,  the  results  at  the  hub  for  other  rotor 
variables  also  differ  from  experimental  data. 

The  exit  Mach  number  and  exit  relative  Mach  number  are  shown  in  Figures 
5.7  and  5.8,  respectively.  Both  curves  follow  the  trends  of  the  experimental  data 
very  closely.  The  absolute  Mach  number  error  at  the  hub  is  due  to  the  near-hub 
deviation  being  under-predicted.  The  maximum  difference  is  -13.83-percent.  The 
relative  Mach  number  is  closer  to  the  experimental  data.  The  maximum  difference 
is  about  4.0-percent  and  occurs  at  the  blade  mid-span. 

The  absolute  and  relative  exit  swirl  velocities  are  discussed  next.  Figure 
5.9  shows  the  exit  swirl  velocity.  This  parameter  is  also  under-predicted  at  the  hub 
because  of  the  near-hub  deviation  error.  The  maximum  difference  is 
-14.12-percent.  The  exit  relative  swirl  velocity,  Figure  5.10,  varies  20.60-percent 
at  the  hub.  Because  this  is  one  of  the  closure  quantities,  the  inaccuracy  in  this 
variable  will  cause  the  entire  flow  field  to  differ  from  the  experimental  data. 

Figures  5.11  and  5.12  show  the  predicted  and  experimental  exit  total  and 
static  pressures,  respectively.  The  total  pressure  is  under-predicted  at  the  hub,  with 
a  maximum  percent  difference  of  —8.84-percent.  At  the  tip,  the  total  pressure  is 
over-predicted.  The  experimental  data  is  much  flatter  than  the  predicted  results. 
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The  exit  static  pressure  is  in  better  agreement  with  experimental  data  with  a 
maximum  difference  of -1.40-percent. 

The  exit  total  temperature,  Figure  5.13,  and  exit  static  temperature,  Figure 
5.14,  each  follow  the  same  trend  as  the  exit  total  and  static  pressures.  The  total 
temperature  is  under-predicted  at  the  hub  by  2.72-percent  and  over-predicted  at  the 
tip.  This  is  a  result  of  the  deviation  being  over-predicted  at  the  hub  and  under¬ 
predicted  at  the  tip.  The  static  temperature  is  predicted  to  within  one-percent  and 
follows  the  experimental  trends  very  closely. 

Figures  5.15  through  5.22  show  the  results  at  the  exit  of  the  stator,  the 
second  blade  row  in  Stage  35.  Because  the  results  of  the  rotor  exit  are  in  error,  the 
inlet  and  exit  flows  of  the  stator  will  also  be  in  error.  Figure  5.15  shows  the 
predicted  relative  total  pressure  loss  coefficient.  The  SLCC  predicts  a  very 
uniform  loss  across  the  span.  The  experimental  results  are  a  factor  of  two  higher 
and  not  nearly  as  uniform.  Unlike  the  rotor,  the  stator  misses  the  exit  flow  angle 
(Figure  5.16)  more  at  the  tip  than  at  the  hub.  This  results  from  the  error  in  the 
flow  deviation  prediction.  The  maximum  error  in  the  exit  flow  angle  is 
-37.75-percent,  or  almost  five  degrees,  at  the  tip. 

Because  the  stator  is  a  stationary  blade  row,  the  absolute  and  relative 
frames  of  reference  are  the  same;  therefore,  no  relative  quantities  are  presented  for 
the  stator.  The  predicted  exit  Mach  number  in  Figure  5.17  tends  to  increase 
toward  the  tip,  as  opposed  to  the  experimental  data  which  decreases  toward  the  tip. 
The  maximum  percent  difference  in  Mach  number  is  almost  31 -percent  at  the  tip. 
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The  stator  is  attempting  to  de-swirl  the  flow  coming  out  of  the  rotor.  This 
can  be  seen  by  comparing  the  stator  exit  swirl  velocity  (Figure  5.18)  with  the  rotor 
exit  swirl  velocity  (Figure  5.9).  The  stator  swirl  velocity  is  approximately  five 
times  lower  than  the  rotor  swirl.  The  predicted  swirl  velocity  is  lower  than  the 
experimental  data.  The  lower  predicted  value  occurred  because  the  deviation  was 
too  low  and  as  a  result  the  flow  was  not  turned  enough. 

Figures  5.19  and  5.20  show  the  stator  exit  total  and  static  pressure, 
respectively.  The  stator  total  pressure  is  in  better  agreement  than  that  of  the  rotor 
exit  total  pressure.  The  maximum  percent  difference  is  only  -6.18-percent.  The 
SLCC  predicts  a  larger  pressure  at  the  tip  than  the  experimental  data.  The  static 
pressure  (Figure  5.20)  is  less  than  the  experimental  data  everywhere  except  at  the 
hub.  The  maximum  percent  error  is -3.30-percent. 

The  stator  exit  total  temperature  and  static  temperature  are  plotted  in 
Figures  5.21  and  5.22,  respectively.  The  total  temperature  trend  is  similar  to  the 
total  pressure,  under-predicted  at  the  hub  and  mid-span  and  over-predicted  at  the 
tip.  The  static  temperature  closely  follows  the  trend  of  the  experimental  data,  only 
slightly  offset.  The  maximum  difference  in  the  total  and  static  temperatures  are 
-2.71-  and  -2.34-percent,  respectively. 

Using  the  averaging  technique  to  convert  radial  distributions  of  flow  field 
variables  into  a  1-D  representation,  the  total  pressure  and  total  temperature 
characteristics  for  Stage  35  were  computed  and  are  presented  in  Figures  5.23  and 


57 


5.24.  The  pressure  characteristics,  VFP,  are  shown  in  Figure  5.23  for  70-,  90-,  and 
100-percent  design  corrected  speeds.  Each  point  of  these  results  is  a  mass  average 
of  a  radial  distribution  of  SLCC  results  at  the  stage  inlet  and  exit.  The  results 
shown  in  Figure  5.23  indicate  that  the  overall  performance  characteristics  follow 
the  experimental  curvature  trend  fairly  well.  The  90-percent  curve  follows  the 
experimental  data  closely,  with  a  maximum  percent  difference  of  only 
4.31-percent.  The  100-precent  speed  line  had  a  maximum  error  of  5.04-percent. 
The  70-percent  speed  line  has  the  highest  maximum  error,  6.42-percent,  at  the 
highest  mass  flow  point.  The  maximum  error  of  6.42-percent  for  the  uncalibrated 
correlations  is  almost  twice  as  accurate  as  the  11.10-percent  difference  between 
experimental  data  and  Stage  35  predictions  with  Rotor  IB  calibration  trends 
included  (illustrated  in  Figure  5.3).  Figure  5.23  clearly  demonstrates  that  even 
though  the  radial  distributions  of  compressor  flow  fields  might  be  predicted  with 
significant  error,  the  mass-averaged  results  of  the  experimental  and  predicted 
results  may  be  quite  acceptable. 

Figure  5.24  displays  the  temperature  characteristics  for  Stage  35.  These 
results  again  show  the  90-percent  speed  prediction  is  better  than  at  the  other  two 
speeds.  The  70-percent  speed  line  has  a  maximum  error  of  1.91-percent  at  the 
highest  mass  flow  point.  The  100-percent  line  has  a  maximum  percent  error  of 
-2.26-percent  at  the  stall  limit.  This  error  is  smaller  than  the  -3.41-percent  error 
shown  in  Figure  5.4  for  the  case  of  Stage  35  with  Rotor  IB  calibration  trends 
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included.  The  percent  error  for  the  mass-averaged  1-D  temperature  characteristic 
representation  (Figure  5.24)  is  again  better  than  the  2-D  representation  of  the  same 
data.  Typically,  the  percent  difference  of  error  in  the  temperature  characteristic  is 
smaller  than  that  of  the  pressure  characteristic. 

The  last  set  of  results  presented  for  Stage  35  are  the  pressure  and 
temperature  characteristics  at  the  hub,  mid-line,  and  tip  for  the  different  speeds. 
The  experimental  data  for  the  tip  is  represented  by  a  diamond,  the  mid-line  a 
circle,  and  the  hub  an  X.  The  results  for  the  tip  are  a  dashed  green  line,  midline  a 
solid  blue  line,  and  the  hub  a  red  dot-dashed  line.  The  pressure  characteristic 
results  are  presented  in  Figure  5.25.  The  predictions  for  all  three  speed  lines  have 
the  tip  with  the  highest  total  pressure  ratio,  the  mid-line  with  the  next  highest 
pressure  ratio,  and  the  hub  with  the  lowest  pressure  ratio.  The  experimental  results 
differ,  however.  At  100-percent  speed,  the  mid-line  has  the  highest  pressure  ratio. 
This  is  because  Stage  35  was  designed  with  a  MCA  blade  profile  for  large  stall 
margin.  Because  the  tip  is  normally  the  first  section  to  stall  on  a  compressor  stage, 
the  tip  was  designed  to  be  less  highly  loaded  at  the  higher  speeds.  This 
characteristic  is  also  present  at  the  90-percent  speed.  At  70-percent  speed,  the  tip 
is  the  most  highly  loaded  section. 

The  maximum  percent  error  for  the  hub  region  is  -6.22-percent  on  the 
100-percent  speed  line.  At  the  midline,  the  70-percent  speed  line  has  the  largest 
error  of  5.59-percent.  The  largest  error  at  the  tip  is  23. 32-percent  on  the 
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90-percent  speed  line.  The  tip  has  the  largest  error  because  the  correlations  used 
were  not  able  to  handle  MCA  blade  profiles  very  well. 

Figure  5.26  shows  the  total  temperature  characteristic  and  has  the  same 
legend  as  Figure  5.25.  The  total  temperature  characteristic  for  Stage  35  exhibits 
different  behavior  than  Rotor  IB.  Stage  35  typically  has  the  highest  temperature 
ratio  at  the  hub  and  the  lowest  ratio  at  the  tip,  whereas,  Rotor  IB  typically  has  the 
opposite  distribution.  The  results  predict  the  tip  to  have  the  highest  temperature 
ratio  and  the  hub  the  lowest  at  all  speeds  for  Stage  35.  This  is  again  due  the 
correlations  which  were  developed  for  DCA  blade  profiles  while  Stage  35  has 
MCA  profiles. 

The  maximum  percent  error  for  the  hub  region  is  -2.82-percent  on  the 
100-percent  speed  line.  In  the  mid-line  region,  the  maximum  error  is 
-3.12-percent  on  the  100-percent  speed  line.  The  100-percent  speed  line  also  has 
the  maximum  percent  error  at  the  tip.  The  error  there  is  6.17-percent.  Typically, 
the  temperature  ratio  percent  error  is  smaller  that  of  the  pressure  ratio  error. 
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6.0  SUMMARY  AND  CONCLUSIONS 


This  chapter  summarizes  the  results  and  conclusions  of  the  work  contained 
in  this  investigation.  The  summary  of  the  results  lists  the  highlights  of  the  results 
of  the  calibration  of  Rotor  IB  and  the  prediction  of  Stage  35  characteristics.  The 
conclusions  discuss  the  lessons  learned  and  the  significance  of  the  results. 

6.1  Summary  of  the  Results 

The  purpose  of  this  thesis  was  to  develop  a  technique  to  generate  stage-by- 
stage,  pre-stall  compressor  performance  characteristics  by  coupling  an  existing 
SLCC  using  open  literature  correlations  for  flow  losses  and  deviations  for  bladed 
regions.  The  method  was  developed,  one  machine  was  calibrated,  and  one 
machine  was  predicted.  The  results  of  applying  the  method  to  the  two  machines 
are  summarized  below: 

•  A  method  to  generate  1-D  and  2-D  compressor  performance 
characteristics  was  developed. 

•  With  calibration  of  relative  total  pressure  loss  coefficient,  flow 
deviation,  and  inlet  and  exit  radial  blade  blockage,  the  entire  flow  field 
can  be  predicted. 

•  A  calibrated  flow  field  for  Rotor  IB  gave  a  1-D  representation  of 
compressor  performance  characteristics  to  within  3-percent. 

•  A  calibrated  flow  field  for  Rotor  IB  gave  a  2-D  representation  of 
compressor  performance  characteristics  to  within  6-percent. 
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•  The  calibration  trends  of  Rotor  IB  imposed  on  Stage  35  gave 
unsatisfactory  results,  therefore  Stage  35  was  predicted  using  the 
original,  uncalibrated  correlations. 

•  An  uncalibrated  flow  field  for  Stage  35  gave  a  1-D  representation  of 
compressor  performance  characteristics  to  within  7-percent. 

•  An  uncalibrated  flow  field  for  Stage  35  gave  a  2-D  representation  of 
compressor  performance  characteristics  to  within  24-percent. 

The  following  section  discusses  the  conclusions  of  what  was  learned  from 
the  results. 

6.2  Conclusions 

Several  important  pieces  of  information  were  learned  from  the  calibration 
and  prediction  of  characteristics  for  Rotor  IB  and  Stage  35.  The  first  conclusion  is 
that,  with  the  four  quantities  used  in  the  calibration,  a  machine  that  does  not  meet 
all  the  constraints  of  the  correlations  can  be  calibrated  if  several  important  flow 
field  quantities  are  known;  in  particular,  radial  distributions  of  loss,  deviation,  total 
pressure,  and  total  temperature. 

Another  piece  of  knowledge  gained  is  that  it  is  easier  to  match  a  1-D 
representation  of  characteristics  than  a  2-D  representation.  For  Stage  35,  the  radial 
distribution  of  the  flow  field  proved  to  be  somewhat  inaccurate  for  several  flow 
field  quantities.  With  the  radial  distribution  mass-averaged  into  a  1-D 
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representation,  the  results  were  in  better  agreement  with  the  mass-averaged 
experimental  1-D  performance  quantities. 

Solidity  is  a  very  important  input  quantity  for  the  SLCC  and  the 
correlations  used  in  this  investigation.  A  look  at  the  Rotor  IB  uncalibrated  results 
shows  that  the  predicted  relative  total  pressure  loss  coefficient  is  significantly  too 
high  for  the  high  solidity  region  at  the  hub.  The  correlations  were  good  for  a 
solidity  of  less  than  two  and  the  hub  of  Rotor  IB  exceeded  this  value.  The 
corrections  applied  to  Rotor  IB,  especially  the  negative  inlet  and  exit  blockage  and 
negative  loss  at  the  hub,  were  able  to  correct  for  the  inaccuracy  of  the  loss  and 
deviation  correlations  at  the  hub. 

Annulus  wall  geometry  also  plays  an  important  part  in  predicting 
characteristics  using  a  SLCC.  The  hub  geometry  for  Rotor  IB  has  a  large  slope 
associated  with  it.  This  large  slope  requires  a  large  curvature  of  the  streamlines 
across  the  blades.  This  leads  to  3-D  flow  effects.  The  correlations  were  not 
designed  to  handle  these  or  any  other  severe  3-D  effects  which  led  to  a  poor 
prediction  of  the  loss  and  deviation  at  the  hub 

The  blade  shape  is  also  identified  as  an  error  contributor  in  the  use  of  the 
deviation  correlations.  The  correlations  used  in  the  SLCC  were  developed  for  the 
double-circular-arc  (DCA)  blade  profile.  Stage  35  blades  use  multiple-circular-arc 
(MCA)  blade  profiles.  The  loss  for  Stage  35  was  predicted  fairly  accurately.  The 
deviation  was  not.  It  is  believed  that  this  blade  shape  difference  is  the  reason  the 
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correlations  could  not  accurately  predict  deviation  for  the  MCA  profiles  of  Stage 
35. 

Another  key  piece  of  knowledge  is  that  the  SLCC  has  much  more  difficulty 
with  convergence  as  the  machine  approaches  the  choked  condition.  The  addition 
of  inlet  and  exit  radial  blockage  also  adds  to  the  difficulty  of  convergence  near 
choke.  A  point  on  the  100-percent  speed  line  of  Rotor  IB  could  not  be  completely 
calibrated  because  any  additional  blockage  would  cause  the  SLCC  solution  to 
diverge. 
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7.0  RECOMMENDATIONS 


Although  the  calibration  of  the  SLCC  and  loss  and  deviation  correlations 
for  Rotor  IB  was  extremely  successful,  the  prediction  of  Stage  35  characteristics 
was  less  accurate.  There  were  also  several  areas  identified  that  could  be  improved 
or  added  to  in  future  work  on  the  same  subject  to  make  the  stage  characteristic 
method  a  more  useful  and  better  predictive  tool.  A  summary  of  recommendations 
for  future  work  is  included  below. 

The  first  recommendation  is  to  use  the  method  to  calibrate  several  more 
machines  that  have  experimental  data  documented.  It  is  apparent  that  the 
calibration  of  the  method  based  on  a  fan  blades  cannot  be  imposed  on  a  core 
compressor  blade  and  get  accurate  results.  Several  fan  rotors  and  stators,  and 
several  core  compressor  rotors  and  stators  should  be  calibrated.  The  calibration 
trends  should  be  investigated  to  see  if  fan  blades  of  similar  families  all  have 
similar  trends  and  core  blades  of  similar  families  have  similar  trends.  A  fan  blade 
or  core  blade  family  should  be  defined  by  blades  with  similar  design  variables 
such  as  stage  loading,  relative  Mach  numbers,  diffusion  factors,  and  solidity.  It 
also  should  be  investigated  to  see  if  rotor  blades  and  stator  blades  have  similar 
calibration  trends. 

The  next  recommendation  is  to  initiate  an  open-literature  investigation  to 
see  if  more  correlations  can  be  found  to  model  blade  rows  with  higher  solidity, 
since  modem  machines  are  tending  toward  higher  solidity. 
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Another  recommendation  is  to  conduct  an  open-literature  investigation  to 
see  if  there  are  correlations  to  handle  additional  types  of  blade  profiles.  The 
correlations  in  the  SLCC  used  in  this  research  were  available  for  NACA  65-(Ai0) 
series  and  double-circular-arc  blade  profiles.  Blading  on  modem  machines  is 
becoming  more  complicated,  and  correlations  should  be  found  with  these  types  of 
profiles  if  this  is  to  truly  be  a  predictive  tool.  If  correlations  cannot  be  found, 
corrections  to  these  correlations  could  be  found  through  calibration  trends  of 
several  machines  with  the  same  types  of  blade  profiles. 

An  alternative  to  the  previous  recommendation  would  be  to  develop 
correlations  for  blades  with  complex  shapes  by  using  3-D  CFD  models.  A 
database  of  performance  variables  for  a  family  of  complex  blades  could  be 
gathered  from  simulations  run  at  on-  and  off-design  conditions.  Correlations  could 
then  be  developed  from  this  database  of  information. 
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Table  4. 1 .  Error  for  Uncalibrated,  Partially  Calibrated,  and  Completely 
Calibrated  SLCC  Data.  (All  Given  in  Percent  Except  for  Loss-Actual  Error) 


Quantity 

Uncalibrated 

Error 

Partially 

Calibrated 

Error 

Calibrated 

Error 

Loss 

-0.15 

0.013 

0.014 

Deviation 

35.56 

-6.48 

-6.58 

Exit  Total 
Pressure 

2.81 

-3.08 

-0.89 

Exit  Static 
Pressure 

3.35 

-2.04 

0.56 

Exit  Total 
Temperature 

-1.36 

-0.68 

0.18 

Exit  Static 
Temperature 

0.70 

-0.43 

0.10 

Exit  Mach 
Number 

-3.87 

-2.14 

-1.93 

Exit  Relative 
Mach  Number 

9.88 

4.20 

3.79 

Ktra  ffiB  t' jSSHBie! 

-11.05 

-8.06 

-5.32 

Exit  Relative 
Swirl  Velocitv 

22.70 

4.78 

3.25 

Table  5.1.  Estimated  Errors  in  Stage  35  Experimental  Data. 


Measurement  Ouantitv 

Estimated  Error 

Mass  Flow 

±  0.66  lbm/s 

Rotational  Sneed 

±  30  rom 

Flow  Anele 

±1.0  dee 

Temperature 

Pi 

o 

r-H 

+1 

Rotor  Inlet  Total  Pressure 

±  0.01  Dsi 

Rotor  Inlet  Static  Pressure 

±  0.04  psi 

Stator  Inlet  Total  Pressure 

±  0.25  osi 

Stator  Inlet  Static  Pressure 

±0.15  Dsi 
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Compressor  Compressor 

Pressure  Ratio  Temperature  Ratio 


Figure  1.1.  Typical  Gas  Turbine  Compression  System.  (Mattingly,  1996) 


Figure  1.2.  Representative  Compressor  Performance  Map. 
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Figure  3.3.  Velocity  Triangle  Nomenclature  Used  in  SLCC. 


Baseline 

“  '  ^Mean  increases 

. co  Terms  increase 

W(width)  increases 


Figure  3.4.  Representative  Loss  Bucket. 


Design  Diffusion  Factor 

Figure  3.5.  P  Factor  for  Minimum  Loss  Calculation 


Figure  3.6  Schematic  of  Lewis,  Hartman,  Miller  Shock  Model 


Figure  3.7.  Schematic  of  Deviation. 


Figure  3.8.  Schematic  of  Streamtube  Geometry  Change  Due 
to  Axial  Velocity  Variation. 


Inlet  Relative  Mach  Number,  M’j 


Figure  3.9.  Variation  of  Deviation  with  Inlet  Relative  Mach  Number  for 
Different  Radial  Locations  on  Double  Circular  Arc  Sections. 
(Robbins,  Jackson,  andLieblein,  1965) 
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Streamlines  Deviation  Angle  Slope 


Solidity,  a 

Figure  3.10.  Deviation  Angle  Slope  Versus  Solidity. 
(Robbins,  Jackson,  and  Lieblein,  1965) 
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2aO 


Figure  4.1.  Schematic  of  NASA  Rotor  IB  Test  Setup.  (Seyler  and  Goselow,  1967) 
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Figure  4.2.  How  Chart  for  SLCC  Calibration  Techniques 
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Figure  4.3.  Radial  Distribution  of  Relative  Total  Pressure  Loss 
Coefficient  including  Experimental  Data,  Uncalibrated  Data,  Partially 
Calibrated  Data  and  Calibrated  Data  for  Rotor  IB. 
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Figure  4.4.  Radial  Distribution  of  Deviation  including  Experimental 
Data,  Uncalibrated  Data,  Partially  Calibrated  Data  and  Calibrated  Data 

for  Rotor  IB. 
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Exit  Relative  Mach  Number  «  Exit  Mach  Number 


♦  Data  —  ■  —Uncalibrated  ——Partial 


Calibrated 


ure  4.5.  Radial  Distribution  of  Exit  Mach  Number  including  Experimental 
Data,  Uncalibrated  Data,  Partially  Calibrated  Data  and  Calibrated  Data 

for  Rotor  IB. 


Figure  4.6.  Radial  Distribution  of  Exit  Relative  Mach  Number  including 
Experimental  Data,  Uncalibrated  Data,  Partially  Calibrated  Data  and 
Calibrated  Data  for  Rotor  IB. 
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Exit  Relative  Swirl  Velocity  (ft/s) 


♦  Data  —  -  —Uncalibrated  ——Partial  - Calibrated 


Figure  4.7.  Radial  Distribution  of  Exit  Swirl  Velocity  including  Experimental  Data, 
Uncalibrated  Data,  Partially  Calibrated  Data  and  Calibrated  Data  for  Rotor  IB. 


Figure  4.8.  Radial  Distribution  of  Exit  Relative  Swirl  Velocity  including 
Experimental  Data,  Uncalibrated  Data,  Partially  Calibrated  Data 
and  Calibrated  Data  for  Rotor  IB. 


87 


Radius  (in.) 

Figure  4.9.  Radial  Distribution  of  Exit  Total  Pressure  including  Experimental 
Data,  Uncalibrated  Data,  Partially  Calibrated  Data  and  Calibrated  Data 

for  Rotor  IB. 


Radius  (in.) 

Figure  4.10.  Radial  Distribution  of  Exit  Static  Pressure  including  Experimental 
Data,  Uncalibrated  Data,  Partially  Calibrated  Data  and  Calibrated  Data 

for  Rotor  IB. 
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Figure  4.11.  Radial  Distribution  of  Exit  Total  Temperature  including 
Experimental  Data,  Uncalibrated  Data,  Partially  Calibrated  Data  and 
Calibrated  Data  for  Rotor  IB. 


Figure  4.12.  Radial  Distribution  of  Exit  Static  Temperature  including 
Experimental  Data,  Uncalibrated  Data,  Partially  Calibrated  Data  and 
Calibrated  Data  for  Rotor  IB. 
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Figure  4.13.  Radial  Distribution  of  Add-Loss  at  100-Percent  Speed 
for  Various  Mass  Flows  for  RotorlB. 
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Figure  4.14.  Radial  Distribution  of  Add-Deviation  at  100-Percent  Speed 
for  Various  Mass  Flows  for  RotorlB. 


90 


Exit  Blockage 


F  221.21 


-4-220.68 


-*-219.38 


Figure  4.15.  Radial  Distribution  of  Inlet  Blockage  at  100-Percent  Speed 
for  Various  Mass  Flows  for  RotorlB. 


Figure  4.16.  Radial  Distribution  of  Exit  Blockage  at  100-Percent  Speed 
for  Various  Mass  Flows  for  RotorlB. 
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Figure  4.17.  Distribution  of  Add-Loss  Averaged  Over  Mass  for  Rotor  IB 
as  a  Function  of  Percent  Span  for  100-,  90-,  70-,  and  50-Percent  Speed. 


Percent  Span 


Figure  4.18.  Distribution  of  Add-Deviation  Averaged  Over  Mass  for  Rotor  IB 
as  a  Function  of  Percent  Span  for  100-,  90-,  70-,  and  50-Percent  Speed. 
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gure  4.19.  Distribution  of  Inlet  Blockage  Averaged  Over  Mass  for  Rotor  IB 
as  a  Function  of  Percent  Span  for  100-,  90-,  70-,  and  50-Percent  Speed. 


Figure  4.20.  Distribution  of  Exit  Blockage  Averaged  Over  Mass  for  Rotor  IB 
as  a  Function  of  Percent  Span  for  100-,  90-,  70-,  and  50-Percent  Speed. 
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4> 


Figure  4.21.  1-D  Pressure  Characteristic  Map  for  Rotor  IB. 

O  Experimental  Data  —  -  Uncalibrated  Results 
—  Calibrated  Results 
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Tip  Results 


Predicted  and  Experimental  2-D  Pressure  Characteristic  Map  for  Rotor  IB 


Tip  Results 


Predicted  and  Experimental  2-D  Temperature  Characteristic  Map  for  Rotor  IB 


Figure  5.3.  1-D  Pressure  Characteristic  Map  for  Stage  35 
with  Rotor  IB  Calibration  Trends  Included. 

O  Data  —  Results  of  Stage  35  with  Rotor  IB  Calibration  Trends 


Figure  5.  4.  1-D  Temperature  Characteristic  Map  for  Stage  35 
with  Rotor  IB  Calibration  Trends  Included. 
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Figure  5.7.  Radial  Distribution  of  Exit  Mach  Number  for 
Stage  35  Rotor. 
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♦  Experimental  Data  —  Uncalibrated  SLCC  Results 


Figure  5.8.  Radial  Distribution  of  Exit  Relative  Mach  Number 
for  Stage  35  Rotor. 
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Figure  5.9.  Radial  Distribution  of  Exit  Swirl  Velocity 
for  Stage  35  Rotor. 


♦  Experimental  Data  —  Uncalibrated  SLCC  Results 


Figure  5.10.  Radial  Distribution  of  Exit  Relative  Swirl  Velocity 

for  Stage  35  Rotor. 
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Figure  5.11.  Radial  Distribution  of  Exit  Total  Pressure 
for  Stage  35  Rotor. 


♦  Experimental  Data  —  Uncalibrated  SLCC  Results 
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Figure  5.12.  Radial  Distribution  of  Exit  Static  Pressure 
for  Stage  35  Rotor. 
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Figure  5.13.  Radial  Distribution  of  Exit  Total  Temperature 
for  Stage  35  Rotor. 


♦  Experimental  Data  —  Uncalibrated  SLCC  Results 


Figure  5.14.  Radial  Distribution  of  Exit  Static  Temperature 
for  Stage  35  Rotor. 
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Figure  5.15.  Radial  Distribution  of  Relative  Total  Pressure  Loss 
Coefficient  for  Stage  35  Stator. 

♦  Experimental  Data  —  Uncalibrated  SLCC  Results 


Figure  5.16.  Radial  Distribution  of  Exit  Flow  Angle  for  Stage  35  Stator. 
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Exit  Swirl  Velocity  (ft/s)  Exit  Mach  Number 


Figure  5.17.  Radial  Distribution  of  Exit  Mach  Number  for 
Stage  35  Stator. 


♦  Experimental  Data  —  Uncalibrated  SLCC  Results 


Figure  5.18.  Radial  Distribution  of  Exit  Swirl  Velocity 
for  Stage  35  Stator. 
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Figure  5.19.  Radial  Distribution  of  Exit  Total  Pressure 
for  Stage  35  Stator. 

♦  Experimental  Data  —  Uncalibrated  SLCC  Results 


Figure  5.20.  Radial  Distribution  of  Exit  Static  Pressure 
for  Stage  35  Stator. 


ITS 


Radius  (in.) 


Figure  5.21 .  Radial  Distribution  of  Exit  Total  Temperature 
for  Stage  35  Stator. 


♦  Experimental  Data  —  Uncalibrated  SLCC  Results 
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Figure  5.22.  Radial  Distribution  of  Exit  Static  Temperature 
for  Stage  35  Stator. 
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Figure  5.23.  1-D  Pressure  Characteristic  Map  for  Stage  35. 
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Figure  5.24.  1-D  Temperature  Characteristic  Map  for  Stage  35. 
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Figure  5.25.  Predicted  and  Experimental  2-D  Pressure  Characteristic  Map  for  Stage  35 
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Appendix  C  -  Summary  of  the  Derivation  of  the 
SLCC  Governing  Equations 

Because  the  SLCC  is  an  integral  part  of  the  method  presented  in  this  thesis, 
it  is  important  to  understand  the  equations  governing  the  code.  This  appendix  will 
describe  the  governing  equations  and  how  the  SLCC  calculates  streamlines 
through  duct  flow  and  through  bladed  regions. 

C.l  Governing  Equations 

The  governing  equations  are  derived  from  the  Navier-Stokes  equations. 
With  the  assumptions  of: 

•  Inviscid 

•  Axisymmetric 

•  No  body  forces 

•  Adiabatic 


•  Steady-state 
the  equations  reduce  to 


r  or  dz 

Continuity 

[C.l] 

dVe  T7  dVe  VrVe  A 
y  0  +y  e  +  '  *  ~o 

dr  z  dz  r 

Circumferential  Momentum 

[C.2] 

v3F-+v3v-=-1^ 

r  dr  z  dz  pdz 

Axial  Momentum 

[C.3] 
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II 

+ 

TO 

** 

1  dP 
p  dr 

Radial  Momentum 

[C.4] 

n 

o 

Energy 

[C.5] 

^  =  0 

Dt 

Entropy 

[C.6] 

P  =  pRT 

Ideal  Gas  Equation  of  State 

[C.7] 

For  use  in  the  SLCC,  the  conservation  equations  were  mapped  to  the  non- 
orthogonal  m  and  /  coordinate  system  shown  in  Figure  3.2.  The  m  direction  is  the 
meridional,  or  streamline,  direction.  The  l  direction  is  the  computing  station 
direction  that  can  be  leaned  to  better  approximate  leading  and  trailing  edges  of 
blades.  Mapping  Equation  [C.2]  to  these  coordinates  gives 

^-K)-o.  [C.8] 

am 

This  states  that  the  angular  momentum  is  constant  along  a  streamline. 

Equation  [C.3]  maps  to 


„  /  2  n\dVm  w  dVm  _ dS  dH  V2cos{(/)  +  y) 

vaw/j)-^=smfr+yy.-^-r-+-+  * 


,  d(tan/?)  cosy 


V2 

— —  tan 2  /? cos  y -  V2  tan 2  ft 
r  dl  sin  (j) 


f 


dS  dH\ 


T - 

dm  dm 


[C.9] 


Equation  [C.4]  maps  to 


dm 


dm  dl  dl 
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Upon  close  inspection  of  Equations  [C.9]  and  [C.10],  they  are  identified  as 
identical  except  for  the  last  terms  with  enthalpy  and  entropy.  By  mapping 
Equations  [C.5]  and  [C.6]  to  the  meridional  direction,  entropy  and  enthalpy  are 
shown  to  remain  constant  along  a  streamline.  Therefore,  the  last  terms  of 
Equations  [C.9]  and  [C.10]  are  identically  zero.  Therefore,  Equations  [C.9]  and 
[C.10]  are  equivalent  to  each  other.  Therefore,  one  of  the  equations  can  be 
discarded.  It  is  decided  to  discard  the  axial  momentum  equation.  The  unknowns 


in  Equation  [C.10]  are  Vm,  /?,  (j),  T,  S,  H,  rc,  and 


The  next  few  steps  will  be 


presented  to  reduce  the  number  of  unknowns  to  the  ones  that  can  be  solved. 

By  using  the  rules  of  differentiation  and  annulus  geometry,  it  is  discovered 

that 

^=(1  tan^tanr)3VM  Kitan<Z)  s™2  ^cos^  rf(tan^)+  sing)  dV^  [Qn] 
dm  K  ^  dz  rc  cosy  dl  cosy  dl 

However,  we  have  introduced  another  unknown,  the  axial  gradient  of  axial 

velocity.  From  continuity  this  can  be  shown  to  be 

9V7  sin  <j)  sin  y  dVm  |  smycos([>Vm  Vm  dp  Vm  sin  </)  cos  ^  sin  0  dVm 
dz  cos(<p  +  y)  dm  cos (<p  +  y)  rc  p  dm  r  cos{(f>  +  y)  dl 


Vm  cos  4  (p  d  (tan  <f>)  Vm  d/ l 
cos  (<p  +  y)  dl  X  dm 


[C.12] 
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This  however  introduces  another  unknown,  — ,  which  by  assuming  perfect  gas 

dm 

and  knowing  that  entropy  and  enthalpy  are  constant  along  a  streamline  can  be 
shown  as 


L^=-m2\  1  dVm  tan2/?sin^ 


p  dm 


V  dm 


[C.13] 


Then,  by  combing  Equations  [C.10],  [C.ll],  [C.12],  and  [C.13],  the  radial 


momentum  becomes 


+a(iK~b(i), 


[C.14] 


where 


i-M: 


^+r)d^l+^k±iliL 

dl  y  dm 


sin  (j)  sin(^  +  /)( 1  +  M%) 


[C.15] 


/  \  2  o  (,  T  )dT.  1  P  dP 

B(l)=  2cos2  c  1 - - - - — — 

W  p  T,  dl  pPt  dl 


[C.16] 


The  continuity  equation  is 
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case  case 

W  =  fdW=  jpv„  cos(0  +  y)MA .  '  [C.  17] 

hub  hub 

With  these  equations  defined,  the  duct  flow  and  bladed  region  solutions  will  be 
described  in  Sections  C.2  and  C.3,  respectively. 

C.2  Duct  Flow 


The  duct  flow  case  involves  no  bladed  region.  The  eight  unknowns  for  this 


case  are:  Vm2,  <fh,  V02, 

Pt2,  Tc  P2,  T2,  pi. 

The  eight  equations  that  solve  these 

unknowns  are: 

^-+A(;K=fl( i) 

at 

Radial  Momentum 

[C.14] 

J~(rVa)=0 

dm 

Circumferential  Momentum 

[C.8] 

case  case 

W  =  [  dW  =  f  pVm  cos(^  +  y)AdA 

Continuity 

[C.17] 

hub  hub 

Dh, 

PD, 

=  0 

Energy 

[C.5] 

^  =  0 

Dt 

Entropy 

[C.6] 

P  =  pRT 

Ideal  Gas  Equation  of  State 

[C.7] 

T=Tt 

+ 

V 

1 

Static  Temperature 

[C.18] 

2cp 

P  =  Pt 

r 

(  t  y-1 

rrt 

Isentropic  Relationship 

[C.19] 
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where  the  radial  momentum,  circumferential  momentum,  continuity,  energy, 
entropy,  and  ideal  gas  equation  of  state  were  defined  in  Section  C.l.  The  three 
closure  parameters  needed  to  resolve  the  flow  field  are  exit  swirl  velocity,  exit 
total  temperature,  and  exit  total  pressure.  By  inspection  of  Equations  [C.5],  [C.6], 
and  [C.8],  it  is  known  that  enthalpy,  entropy,  and  angular  momentum  are 
conserved  along  a  streamline.  Therefore,  the  exit  swirl  velocity  can  be  resolved  by 
integrating  the  circumferential  momentum  equation  (Equation  [C.8])  to  give 


fr.\ 

-  vn. 

yr2  J 


[C.20] 


Then,  the  exit  total  temperature  is  found  by  integrating  the  energy  equation 
(Equation  [C.5])  to  give 


T  =T 

Lt  2  xtl' 


[C.21] 


Finally,  by  integrating  the  entropy  equation  (Equation  [C.6]),  the  exit  total 
pressure  relation  becomes 

p*=p«-  rc-22J 

With  these  closure  relations  and  the  conservation  equations  listed  at  the  beginning 
of  this  section,  the  entire  flow  field  for  a  non-bladed  region  can  be  resolved. 


C.3  Bladed  Region 


Resolving  the  flow  field  in  the  bladed  region  is  more  complicated.  The 
unknowns  and  several  of  the  governing  equations  are  the  same  as  the  duct  flow 
case.  The  closure  relations  are  the  only  equations  that  change.  Therefore,  the 
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radial  momentum  (Equation  [C.14]),  continuity  equation  (Equation  [C.17]),  ideal 
gas  equation  of  state  (Equation  [C.7]),  static  temperature  equation  (Equation 
[C.18]),  and  the  isentropic  relationship  (Equation  [C.19])  are  all  the  same  for  the 
duct  flow  and  bladed  region  cases.  The  new  closure  relations  are  developed 
because  entropy,  enthalpy,  and  angular  momentum  are  no  longer  constant  along  a 
streamline  and  must  be  specified.  Because  of  the  complex  flow  involved  when 
crossing  a  bladed  region,  the  SLCC  does  not  directly  model  the  bladed  regions. 
For  the  closure  relations  to  be  specified,  loss  and  deviation  across  a  bladed  region 
must  be  known.  For  the  SLCC  used  in  this  investigation,  the  loss  and  deviation 
were  specified  using  loss  and  deviation  correlations.  Loss  and  deviation  are 
defined  in  Section  3.1. 

With  deviation  specified  by  the  correlations  and  the  velocity  triangles  and 
exit  streamline  angles  known  at  the  exit  radius,  the  first  of  the  closure  relations, 
exit  swirl  velocity,  can  be  calculated  from 


tan  = 


[C.23] 


cos  = 


[C.24] 


K2=Vm2  tan  /?2  cos  02  +  ^2  • 


[C.25] 


The  second  closure  relation  uses  two  relationships  for  work  to  calculate  the 
exit  total  temperature.  By  equating  the  Euler  turbomachinery  equation, 
w  =  U2Ve2-UlV6l,  [C.26] 
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to  the  first  law  of  thermodynamics  with  the  assumptions  of  steady-state  and 
isentropic  flow  with  no  heat  addition, 

w  =  ht2  -  hn ,  [C.27] 

the  equation  for  exit  total  temperature  is  determined  to  be 


Tt2=Tn+U2 


V  -—V 
v  62  v ' 


6 1 


1 


[C.28] 


Jcp 


The  final  closure  relation  is  used  to  calculate  the  exit  total  pressure.  This 
relationship  uses  the  thermodynamic  quantity  rothalpy,  defined  as 

/  =  A  +  ^1-  — .  [C.29] 

2  2 

For  rothalpy  to  be  conserved  in  a  moving  passage,  the  flow  must  be  steady  in  the 
rotating  frame,  no  work  can  be  performed  in  the  rotating  frame,  and  no  heat 
transfer  can  occur  to  or  from  the  flow.  Then,  using  the  isentropic  relations,  the 
relative  total  pressure  loss  coefficient,  and  the  conservation  of  rothalpy,  the  exit 
total  pressure  is  determined  to  be 


^°I2  —  ^tl 


where 


f  n»  A 

r  t\ 

P* 

V  '2  M 


( T 
1ti 

\Tn  j 


21 

7-1 


1- 


(  p<  \ 

r  n 
P' 

V  <2  Jid 


w' 


1- 


1+ 


r- 1 


_r_  \ 
r-i 


[C.30] 


\  +  — — -M 


n 


\ 


-1 


7 

r- 1 


[C.31] 
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and 


M  =  U'  -  [C.32] 

" 

With  the  three  closure  equations  defined,  Equations  [C.25],  [C.28],  and  [C.30],  the 
conservation  equation  listed  at  the  top  of  this  section  can  be  solved  to  develop  the 
flow  field  for  a  bladed  region. 
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Compressor  Performance  Characteristic  Forms 
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Appendix  D  -  Compressor  Performance  Characteristic 

Forms 

Compressor  performance  characteristics  can  be  presented  in  different 
forms  depending  on  user  requirements  and  application.  The  following  will 
provide  a  discussion  of  four  different  forms  available  in  the  current  investigation. 
These  fomrs  are  all  discussed  in  more  detail  by  Garrard  (1995). 

The  first  form  is  considered  the  classical  form.  The  stage  flow  coefficient 
is  defined  as 

<b  =  -  [D.l] 

U 

and  the  pressure  and  temperature  coefficients,  respectively,  are  defined  as 
¥ p=PR  IP-2] 

'Fr=77?-1.  [D-3] 


The  second  form  is  defined  by  using  the  concepts  of  Mach  number,  flow 
function  and  critical  reference  state  and  are  given  by 


.  wUiv/O 

[D.4] 

V  =prt(nrcJ 

[D.5] 

vT  =(TR-iXmcorr)2. 

[D.6] 

where 
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